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Abstract 


Infrared  Multiphoton  Absorption  (IRMPA)  was  used  to  produce 
populations  of  vibrationally  excited  1,1,2-trifluoroethane,  which  were 
characterized  by  optoacoustlc  measurements  of  absorbed  laser  power 
and  collision  free  decomposition  yields.  The  measurements  were 
accurately  fitted  with  a  Master  Equation  that  included  Quack's 
theory  of  IRMPA,  three  RRKM  unimolecuiar  reaction  channels,  and 
collisional  energy  transfer.  The  highly  constrained  adjustable 
parameters  indicate  that  the  optical  coupling  matrix  elements  are 
dramatically  reduced  in  magnitude  near  reaction  threshold  energies 
where  vibrational  anharmonicity  becomes  important.  Observed 
infrared  fluorescence  from  the  excited  molecules  is  in  excellent 
agreement  with  Master  Equation  predictions  and  it  was  used  to 
monitor  collisional  deactivation  of  the  excited  molecules.  The  energy 
transfer  exhibits  a  weak  vibrational  energy  dependence  and  no 
detectable  temperature  dependence.  Experiments  were  undertaken 
to  ascertain  the  effect  of  vibrational  energy  on  bimolecular  reactions 
of  TFE,  but  no  such  effects  have  yet  been  observed. 

A  Monte  Carlo  method  is  described  for  efficient  multi¬ 
dimensional  integration  not  restricted  to  hyper -dimensional 
rectangles,  but  applied  to  more  complicated  domains.  When  known, 
the  boundaries  of  an  arbitrary  integration  region  can  be  used  to 
define  the  sampling  domain,  resulting  in  sampling  with  unit 
efficiencyiFor  sums  of  states,  the  method  can  include  all  molecular 
degrees  of  meedom,  coupled  in  any  fashion.  The  method  is  highly 
efficient  and  the  results  show  that  the  off-diagonal  anharmonicities 
produce  significant  increases  in  the  calculated  sum  of  vibrational 
states. 
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Appendix  A:  Vibrationally  excited  populations  from  IR-multiphoton 
absorption.  I.  Absorbed  energy  and  reaction  yield 
measurements 

Appendix  B:  Vibrationally  excited  populations  from  IR-multiphoton 
absorption.  II.  Infrared  fluorescence  measurements 

Appendix  C:  Vibrationally  Excited  Populations  from  IR-Multiphoton 
Absorption.  3.  Energy  Transfer  between  1,1,2-Trifluoroethane 
and  Argon 

Appendix  D:  Sums  of  Quantum  States  for  Nonseparable  Degrees  of 
Freedom.  Multidimensional  Monte  Carlo  Integration 


Introduction 


Energetic  materials  that  combine  high  performance  with  low 
sensitivity  have  long  been  sought.  Unfortunately,  many  problems 
are  associated  with  such  materials,  including  irregular  burning, 
premature  detonation,  and  toxicity  Accidental  initiation  of  materials 
with  inherent  molecular  sensitivity  has  been  one  of  the  most 
persistent  and  formidable  problems.  This  serious  problem  motivates 
much  research  on  how  energetic  material  sensitivity  and 
performance  are  related  to  the  molecular  properties  of  the 
ingredients.1  Despite  much  work  on  the  relationships  connecting 
performance  and  sensitivity  to  molecular  properties,  the  detailed 
reaction  sequences  leading  to  detonation  are  not  fully  understood  2 
Other  factors,  such  as  aging  and  environmental  conditioning  (even 
hydrogen  bonding!3),  further  complicate  the  picture,  because  they  can 
significantly  alter  the  sensitivity  of  energetic  materials 

The  focus  of  the  ARO-supported  work  carried  out  in  this 
laboratory  has  been  to  investigate  the  fundamental  properties  of 
highly  vibrationally  excited  molecules  (HVEMs)  related  to  energetic 
materials  Vibrational  excitation  is  produced  in  all  high  temperature 
systems,  including  the  extremely  hostile  environments  associated 
with  energetic  material  combustion,  detonation,  and  deflagration. 
Little  is  known  about  the  energy  transfer  properties,  spectroscopy, 
and  chemical  reactions  of  highly  vibrationally  excited  molecules.  It  is 
possible  that  HVEMs  may  exhibit  enhanced  bimolecular  reaction  rates 
and  unusual  chain  branching  reactions. 

In  this  project,  we  have  used  a  prototype  molecule  to 
investigate  methods  for  preparing  HVEMs  and  characterizing  their 
population  distributions.  This  phase  of  work  has  been  very  successful 
and  we  have  completed  perhaps  the  most  detailed  analysis  ever 
carried  out  on  preparation  of  HVEMs  by  infrared  multiphoton 
absorption  (IRMPA).  In  addition,  we  used  the  HVEMs  prepared  by 


1RMPA  to  determine  the  vibrational  energy  transfer  properties  of  the 
HVEMs.  At  present,  we  are  engaged  in  a  search  (so  far  unsuccessful) 
for  vibrational  enhancement  of  bimolecular  reactions. 

Infrared  multiphoton  absorption  (IRMPA)  is  a  powerful  method 
for  producing  large  populations  of  HVEMs  through  the  use  of  a  high 
power  CO2  laser  During  the  excitation  process,  some  molecules 
absorb  many  photons,  while  other  molecules  only  absorb  a  few, 
resulting  in  broad  population  distributions  controlled  by  the  laser 
fluence  The  first  necessity  was  to  develop  methods  for  determining 
the  population  distributions 


Population  Distribution 

To  determine  the  population  distribution,  we  first  selected  a 
particularly  favorable  prototype  molecule  (1,1,2-trifluoroethane,  le. 
"TFE")  and  used  a  combination  of  experimental  techniques  in 
conjunction  with  a  detailed  Master  Equation  model,  as  summarized 
in  Appendix  A.  The  basic  idea  was  to  take  several  experimental 
measures  of  the  population  distribution,  fit  the  observations  with  a 
single  theoretical  model,  and  then  use  the  model  to  fill  in  the  gaps  in 
our  knowledge 

The  first  experimental  measure  was  the  average  amount  of 
energy  absorbed  from  tne  laser;  this  gives  the  average  energy  of  the 
population  distribution.  The  second  experimental  measure  was  the 
decomposition  yield  as  a  function  of  laser  fluence.  Because  only  the 
most  chemically  energetic  molecules  decompose,  this  was  a  measure 
of  the  high  energy  tail  of  the  distribution.  The  Master  Equation 
model  included  collisions,  optical  pumping,  stimulated  emission,  IRF, 
three  RRKM  unimolecuiar  reaction  channels,  etc.  and  it  was  possible 
to  fit  the  data  with  relatively  little  empiricism,  as  described  in 
Appendix  A. 

To  test  the  accuracy  of  the  population  distribution  as  described 
wyh  the  Master  Equation,  we  carried  out  infrared  emission 


experiments  in  which  we  determined  the  IRF  intensity  from  the  c-H 
modes  as  a  function  of  laser  fluence,  which  determines  the 
population  distributions.  As  shown  in  Appendix  B,  a  comparison 
between  the  experiments  and  IRF  Intensities  predicted  by  the  Master 
Equation  model  showed  excellent  agreement,  demonstrating  that  the 
combination  of  techniques  and  the  Master  Equation  model  gave 
accurate  population  distributions. 


Energy  Transfer 

The  first  application  of  the  population  distributions  was  the 
determination  of  the  energy  transfer  properties  of  TFE,  the  prototype 
molecule,  in  collisions  with  argon,  as  described  in  Appendix  C.  In 
these  experiments.  IRF  from  TFE  was  monitored  as  a  function  of 
time.  The  intensity  was  observed  to  decay  with  time  constants  that 
depended  on  laser  fluence  and  on  the  gas  mixture  composition.  High 
laser  fluences  produce  initial  populations  with  higher  excitation,  and 
the  final  temperature  is  governed  by  the  relative  amount  of  TFE  in 
^he  gas  mixture.  Thus,  by  varying  the  gas  composition  and  laser 
fluence  independently,  we  were  able  to  determine  the  dependence  of 
<AE>d  (the  average  energy  lost  in  deactivating  collisions)  on  both 
excitation  energy  and  bath  temperature.  This  was  the  first  such 
measurement  using  TFE  and  the  results  appear  to  be  quite 
reasonable,  when  compared  with  similar  molecules. 

A  continuation  of  the  energy  transfer  experiments  is  currently 
underway.  Our  objective  here  is  to  compare  the  time-resolved 
optoacoustic  technique  used  by  Gordon  and  coworkers  to  the  IRF 
method  mentioned  above  This  method  is  related  to  time-dependent 
thermal  lenstng  and  is  one  of  several  photothermal  techniques  that 
can  be  applied  to  energy  transfer.  In  essence,  the  energy  transfer 
produces  rapid  heating  from  the  deposited  laser  energy  The  heating 
produces  a  compression  wave  that  travels  outward  from  the  laser 
beam.  Associated  with  the  compression  wave  is  an  expansion  wave 
with  amplitude  dependent  on  the  rate  of  V-T  energy  transfer;  thus 


the  ratio  of  the  wave  amplitudes  is  related  to  the  time  constant  for 
energy  transfer.  Either  the  theory  of  Bailey  et  al.,  or  the  more 
complete  theory  of  Barker  and  Rothem  can  be  used  to  describe  the 
pressure  waves  In  this  energy  transfer  system 


Bimolecular  Reactions 

We  have  expended  a  major  effort  over  the  past  year  to  detect 
and  measure  changes  In  bimolecular  reactivity  of  TFE  due  to 
vibrational  excitation.  In  these  experiments,  we  have  attempted  to 
observe  the  reactions  of  H-atoms  and  O-atoms  with  IRMPA-exclted 
TFE.  So  far,  we  have  observed  no  enhancement,  but  we  plan  to 
make  several  more  attempts,  before  concluding  the  experiments. 

Before  beginning  the  experiments  on  bimolecular  reactions,  we 
first  tested  our  new  laboratory  laser  (Lumonics  Hyperex  400)  to 
determine  whether  it  produces  populations  similar  to  those  produced 
by  the  Tachisto  and  Lumonics  103  lasers  originally  used  by  us  at  SRI. 
New  optoacoustic  measurements  were  within  10?  of  the  original  data, 
and  decomposition  yield  measurements  were  also  reproduced  within 
the  experimental  uncertainties  (3?  yield,  compared  with  4?  yield 
obtained  earlier  at  0.9  J  cm-2).  These  levels  of  agreement  are 
surprisingly  good,  because  the  new  measurements  were  carried  out 
with  a  different  laser,  new  power  meter,  different  optoacoustic 
microphone,  and  new  quadrupole  mass  spectrometer. 

The  reaction  that  we  are  currently  Investigating  is  that  of  H- 
atoms  with  TFE  to  produce  H2  and  HF 

H  ♦  CF2H-CFH2  -  H2  «■  C2F3H2  (3a) 

-  HF  *  C2F2H3  (8b) 

We  estimate  that  both  reaction  channels  have  A-factors  of  about 
10'10  cm3  s’1  and  activation  energies  of  about  340  to  400  k  J  mol"1 


and  we  are  in  the  process  of  measuring  the  thermal  rate  constants 
to  confirm  these  estimates. 

If  vibrational  energy  is  important  in  the  reaction,  it  may 
increase  (or  decrease!)  the  total  reactivity  of  TFE  and/or  modify  the 
branching  ratio  for  the  two  channels.  H-atoms  (or  O-atoms)  are 
produced  with  a  microwave  discharge  in  a  flowing  argon  or  helium 
carrier  gas;  the  atom  concentration  is  determined  by  titration  with 
NO2  (the  end-point  is  determined  mass  spectrometrically): 

H  ♦  NO2  -*  OH  ♦  NO  (10) 

This  reaction  is  followed  by  several  subsequent  reactions,  giving  a 
stiochiometry  factor  somewhat  larger  than  unity. 

The  actual  H-atom  (or  0-atom)  concentration  is  needed  only 
when  using  H-atom-rich  conditions,  as  we  have  attempted  to  do,  thus 
far  Unfortunately,  we  have  not  been  able  to  produce  sufficient 
concentration  of  H-atoms  under  the  low  pressure  conditions  needed  to 
avoid  collisional  deactivation  of  the  excited  molecules.  Under  these 
conditions,  we  also  measured  the  reaction  rate  of  H-atoms  with 
ethylene  to  verify  the  H-atom  concentrations.  This  experiment 
showed  an  even  lower  [H],  indicating  that  NO2  was  probably  also 
being  lost  due  to  photolysis  from  the  microwave  discharge. 

To  avoid  these  problems,  we  must  use  TFE-rich  conditions, 
where  the  concentration  of  vibrattonally  excited  TFE  is  in  great 
excess  over  [H].  These  conditions  are  also  desirable  for  another 
reason.  H-atoms  are  likely  to  deactivate  rapidly  the  excited  TFE,  and 
the  low  [H]  will  reduce  the  importance  of  this  process  (we  plan  to 
measure  this  deactivation  rate).  A  major  advantage  of  Infrared 
multiphoton  absorption  is  that  it  is  relatively  easy  to  produce  large 
concentrations  of  vibrationally  excited  molecules.  By  using  30  mTorr 
of  TFE,  we  will  produce  =  1015  cm'3  excited  TFE  molecules,  much 
larger  than  the  =  1013  cm-3  H-atoms  we  are  currently  producing.  To 
work  with  low  [H],  we  will  use  resonance  absorption  and/or 
resonance  fluorescence  at  the  lyman-a  line  in  future  experiments. 


Despite  the  low  H-atom  concentrations,  we  have  carried  out 
experiments  by  using  mass  spectrometry  to  monitor  [TFE]  and  the 
production  of  reaction  products.  Thus  far,  we  have  not  observed  any 
significant  variation  due  to  blmolecular  reaction  in  the  HF  production, 
or  in  loss  of  TFE  when  the  molecule  is  irradiated  at  moderate 
fluences.  These  results  have  been  disappointing,  but  not  exceptional, 
since  our  conditions  have  not  been  optimal.  Moreover,  some 
polyatomic  reaction  systems  have  shown  a  dependence  on  vibrational 
energy,  but  others  have  not.4  Our  experiments  are  still  underway, 
and  monitoring  H-atom  concentrations  in  the  presence  of  excess 
[TFE(vib)]  is  likely  to  be  a  more  sensitive  way  of  observing  changes  in 
the  reaction  rate.  Regardless  of  whether  an  effect  is  ultimately 
observed,  such  experiments  could  not  even  be  attempted  without 
IRMPA  production  of  excited  molecules. 


Theoretical  Work 

As  part  of  the  analysis  of  the  population  distributions,  a  full 
Master  Equation  treatment  was  needed.  Although  much  of  the 
computer  code  had  already  been  written,  we  found  it  necessary  to 
include  M.  Quack's  statistical-dynamical  theory  of  infrared 
multiphoton  excitation,  as  described  in  Appendix  A.  Although  most 
of  the  implementation  of  this  theory  was  straightforward,  we 
eliminated  many  of  the  approximations  that  were  originally  used  by 
Quack.  For  example,  we  used  exact  counts  of  vibrational  states, 
rather  than  the  semiclassical  approximations  used  by  the  original 
authors. 

In  addition  to  the  implementation  of  Quack's  theory,  we  also 
have  carried  out  original  research  on  the  sums  and  densities  of 
states  of  HVEMs  (Appendix  D).  In  real  molecules,  all  degrees  of 
freedom  (DOF)  are  coupled,  as  distinguished  from  the  usual  separable 
harmonic  and  anharmonic  models.  For  separable  DOF,  the  sums  and 
densities  of  states  are  easily  counted  exactly,  through  use  of  the 
Stein -Rabinovitch  algorithm.  However,  when  molecules  are  highly 


excited,  separability  breaks  down  and  coupling  among  the  DOF 
becomes  dominant.  Our  research  on  molecules  with  coupled  DOF  was 
motivated  by  the  need  to  determine  exact  densities  of  states  for 
HVEMs 

To  treat  coupling  among  all  DOF  is  a  major  challenge,  because 
the  couplings  include  the  usual  diagonal  and  off-diagonal  vibrational 
anharmonicities,  couplings  between  vibrations  and  rotations,  corlolis 
effects,  Fermi  resonances,  etc.  For  different  theoretical  applications 
different  ones  of  these  couplings  may  be  of  interest.  Our  approach 
was  to  develop  a  Monte  Carlo  method  for  treating  these  systems, 
because  such  methods  are,  in  principle,  completely  general. 

The  Monte  Carlo  method  we  developed  Is,  we  believe,  a  new 
approach  to  multi-dimensional  integration.  We  employed  it  in 
discrete  state-space  for  the  purpose  of  counting  sums  of  states,  but  it 
can  also  be  used  with  continuous  variables.  The  details  are  given  in 
Appendix  D. 


To  demonstrate  the  calculation  of  sums  of  states  for  non- 
separable  DOF,  we  used  the  spectroscopic  parameters  for  H2O  and  for 
CH2O  to  calculate  the  vibrational  sums  of  states  We  used  the  actual 
spectroscopic  data,  not  a  theoretical  model1  The  results  showed  that 
the  Monte  Carlo  method  we  developed  is  efficient  and  practical.  The 
resulting  sums  of  states  showed  that  separable  models  for  the  same 
molecules  underestimate  the  actual  sums  of  states  by  significant 
factors,  which  could  affect  theoretical  interpretations  in  some  cases. 

Future  work  on  development  of  this  method  will  treat 
couplings  between  vibrations  and  rotations,  which  are  crucially 
important  in  unimolecular  reactions  near  threshold 
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Conclusions 


The  completed  phases  of  this  research  were  highly  successful. 
The  characterization  of  the  population  distribution  of  excited  TFE 
molecules  was  probably  the  most  complete  ever  carried  out  on  any 
molecule  other  than  SF6  The  IRF  emission  test  of  the  distribution 
function  showed  excellent  agreement  with  theory,  and  the  IRF  was 
then  used  to  measure  energy  transfer  between  excited  TFE  and 
argon  collider  gas. 

The  difficult  experiments  to  determine  the  effect  of  vibrational 
energy  on  blmolecular  reactions  are  still  underway  and  no  definitive 
conclusion  has  been  reached,  as  yet.  Using  mass  spectrometry  to 
monitor  TFE  concentrations  and  reaction  products,  we  have  not 
observed  any  influence  of  vibrational  energy.  However,  conditions 
are  not  optimum,  since  the  H-atom  concentrations  are  too  low  to 
provide  a  good  test.  We  will  use  Lyman-a  photometry  and/or 
resonance  fluorescence  to  monitor  H-atom  concentration  in  future 
experiments.  Under  these  conditions,  effects  are  more  likely  to  be 
observed. 

The  theoretical  efforts  have  been  very  successful,  leading  to 
improvements  in  the  collislonal/photophystcal  Master  Equation  used 
in  analysis  of  experimental  results,  as  well  as  in  development  of  a 
new  Monte  Carlo  multidimensional  integration  scheme  and 
evaluation  of  sums  of  states  for  coupled  vibrational  degrees  of 
freedom. 
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The  molecule  1,1,2-trifluroethane  (TFE)  was  used  in  experiments  to  determine  the  population 
distribution  of  excited  molecules  produced  by  infrared  multiphoton  absorption  induced  by  high 
power  TEA  C02  lasers  operating  at  1079.85  cm-1  [9  6  ^m  R  (22)  line].  Optoacoustic 
measurements  of  absorbed  laser  power  provided  a  measure  of  the  mean  energy  of  the  population 
distribution,  while  very  low  pressure  photolysis  measurements  of  the  collision-free  decomposition 
yield  gave  information  about  the  high-energy  tail  of  the  distribution.  The  experimental  results 
were  accurately  simulated  using  a  Master  Equation  model  that  incorporated  Quack’s  statistical- 
dynamical  theory  of  infrared  multiphoton  absorption  (cases  B  and  C),  RRKM  unimolecular 
reactions  (three  channels),  and  collisional  energy  transfer.  The  computer  simulations  included 
known  TFE  molecular  properties  and  only  four  adjustable  parameters,  which  were  very  highly 
constrained  in  order  to  fit  the  experimental  data.  From  the  simulations,  we  conclude  that  the 
optical  coupling  matrix  elements  are  dramatically  reduced  in  magnitude  for  energies  above  the 
reaction  thresholds.  This  effect  is  symptomatic  of  the  vibrational  anharmonicity  due  to  the 
presence  of  the  reaction  channels,  even  in  molecules  that  have  not  yet  reacted,  resulting  in 
vibrational  frequency  shifts  of  the  absorption  lines  out  of  resonance  with  the  laser  line.  This  effect 
is  expected  to  be  present  and  observable  in  other  highly  vibrationally  excited  molecules. 


INTRODUCTION 

A  little  more  than  ten  years  ago.  great  interest  was 
aroused  by  the  first  observations  of  chemical  reactions  in¬ 
duced  by  infrared  lasers.1  The  observations  could  only  be 
explained  if  isolated  molecules  absorbed  many  infrared  pho¬ 
tons.  leading  to  bond  fission  and  molecular  elimination  reac¬ 
tions.  Because  not  all  irradiated  molecules  absorb  the  same 
number  of  photons,  the  population  distributions  are  pro¬ 
duced  with  substantial  energy  dispersions. 

Desptte  the  broad  population  distributions,  IR-muIti- 
photon  absorption  has  great  potential  as  a  useful  laboratory 
tool  because  of  two  features:  ( 1 )  many  molecules  absorb  light 
from  powerful  infrared  (e  g.,  CO:j  lasers  and  therefore  the 
effect  is  general,  and  (2)  almost  all  irradiated  molecules  are 
excited  to  relatively  high  vibrational  energies.  By  varying 
laser  power,  wavelength,  and  other  parameters,  the  average 
energyof  the  population  can  be  controlled,  much  as  the  aver¬ 
age  energy  can  be  controlled  in  thermal  systems  by  changing 
the  temperature.  Until  the  population  distributions  are  bet¬ 
ter  characterized,  however,  the  usefulness  of  IRMPA  is  li¬ 
mited,  despite  its  potential  as  a  laboratory  tool. 

Although  hundreds  of  publications  have  described 
IRMPA  and  several  excellent  reviews  have  been  published,1 


"  All  of  the  experimental  work  waa  carried  out  at  SRI  International  and  was 
supported  by  the  U  S  Army  Research  Office. 
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the  details  of  the  population  distributions  in  large  molecules 
are  still  not  well  known,  but  only  have  been  inferred  from 
limited  measurements2  '’  and  from  theoretical  calculations. 
In  this  paper,  two  experimental  techniques  (VLP<I>  measure¬ 
ment  of  decomposition  yields4  and  optoacoustic  measure¬ 
ment  of  absorbed  laser  poweri  were  employed  along  with  a 
Master  Equation  theoretical  model6  to  obtain  improved  esti¬ 
mates  of  the  population  distributions  produced  by  IRMPA. 
By  using  more  than  one  experimental  technique,  the  theo¬ 
retical  modeling  is  highly  constrained,  leading  to  reliable 
estimates  for  the  population  distributions.  A  third  expen- 
mental  technique,  infrared  fluorescence  spectroscopy,  was 
also  applied  in  this  study  and  it  shows  excellent  consistency 
with  the  other  techniques,  but  it  has  not  been  independently 
calibrated  and  so  it  is  descnbed  separately  in  the  following 
paper.7 

1.1.2-tnfluorethane  iTFE)  was  chosen  for  the  expen- 
ments  for  several  reasons: 

1 1 1  A  vibrational  assignment"  has  been  made  and  TFE 
has  a  strongly  absorbing  C-F  stretch  mode  near  1076 cm'1, 
suitable  for  pumping  with  a  CO,  laser. 

(2)  TFE  has  simple  unimolecular  reactions  (HF  molecu¬ 
lar  eliminations  by  three  reaction  pathways),  which  have  al¬ 
ready  been  thoroughly  investigated  in  thermal  and  chemical 
activation  experiments.’  These  molecular  elimination  reac¬ 
tions  do  not  lead  to  secondary  chemistry,  simplifying  inter¬ 
pretation  of  the  results. 
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CF:H-CFH,— CF,  =  CH,  HF  iA| 

—cu-  or  trans-CFH  =  CFH  -t-  HF  |B) 
-*:CF-CFH2  +  HF 

— <is-  or  tranr-CFH  =  CFH  -t-  HF,  (C) 

1 3)  The  internal  and  external  rotations  in  TFE  are  rela¬ 
tively  massive  and  they  have  low  symmetry,  producing  a 
high  density  of  optically  accessible  vibrational  states,  which 
enhances  the  efficiency  of  IRMPA-  Moreover,  the  wide  dis¬ 
persion  of  the  vibrational  frequencies  produces  a  “filled-tn” 
density  of  states,  with  few  regions  of  low  state  density,  except 
at  very  low  energy. 

(4|  The  C-H  stretch  mode  frequencies  are  conveniently 
separated  from  the  frequencies  of  other  modes  in  TFE,  re¬ 
ducing  the  spectroscopic  resolution  requirements  in  the  in¬ 
frared  fluorescence  experiments.’ 

1 5 1  The  bimolecular  reactions  of  TFE  with  atoms  and 
simple  free  radicals  are  expected  to  be  straightforward  and 
suitable  for  studies  of  the  effects  of  vibrational  energy  on 
bimolecular  reactivity. 

The  usefulness  of  TFE  has  been  borne  out  in  the  experi¬ 
ments  reported  here.  High  IRMPD  yields  were  obtained  and 
the  experimental  observations  were  entirely  in  accord  with 
the  qualitative  expectations.  The  results  are  highly  consis¬ 
tent  and  they  lead  to  better  knowledge  of  the  properties  of 
highly  vibrationally  excited  molecules,  as  described  in  the 
following  sections. 

EXPERIMENTAL 

The  VLP  d>  technique  for  measurement  of  collision-free 
IR  multiphoton  decomposition  ( IRMPD)  yields  has  been  de¬ 
scribed  elsewhere  in  detail,4  except  for  a  few  modifications. 
Gases  at  low  pressure  ( <  20  mTorr)  are  irradiated  in  a  2  cm 
diam  by  10  cm  long  Knudsen  cell  equipped  with  KC1  win¬ 
dows.  The  gas  escapes  through  an  aperture  at  a  rate  depen¬ 
dent  only  on  the  mean  molecular  velocity,  the  cell  volume, 
and  the  area  of  the  aperture.  In  the  present  experiments,  the 
escape  aperture  was  fabricated  from  a  Teflon-bore  greaseless 
stopcock  in  order  to  have  continuous  control  of  the  residence 
time.  Since  the  stopcock  could  not  be  reproducibly  reset  to 
known  aperture  areas,  1  Torr  Baratron  capacitance  manom¬ 
eter  (model  223A)  was  used  to  monitor  gas  pressure  in  the 
cell  directly.  The  residence  time  r  is  given  by10 

r  =  PVuaS0/(RTFA  (,  ill 

where  Pis  the  observed  pressure,  K€„  is  the  cell  volume,  .V', 
is  Avogadro’s  number,  R  is  the  gas  law  constant.  T  the  tem¬ 
perature,  and  Fa  is  the  measured  gas  flow  rate  i  measured  by 
timing  the  rate  of  pressure  increase  in  a  calibrated  volume!; 
the  escape  rate  constant  is  r  ~  *. 

The  composition  of  the  gas  escaping  from  the  cell  is 
monitored  by  a  modulated  molecular  beam  mass  spectrom¬ 
eter,  as  described  elsewhere."  For  this  purpose,  a  Finnegan 
S  pec  t  rase  an  400  quad  ru  pole  mass  spectrometer  was  used 
with  an  SRI-built  “peak  picker"  for  selection  of  specific 
mass  peaks.  The  output  of  the  quadrupole  was  amplified 
with  a  PAR  model  1 8 1  charge-sensitive  preamplifier  and  the 
modulated  signal  was  isolated  from  background  with  an  Ith- 


aco  model  391  lock-in  amplifier,  locked  to  the  200  Ha  tun¬ 
ing-fork  chopper  (Bulova  Time  Products).  For  father  signal 
averaging,  the  output  from  the  lock-in  was  captured  and 
stored  by  a  Nicolet  model  1072  signal  averager,  which  was 
interfaced  with  a  DEC  LSI  1 1/2  computer  for  storage  of  the 
data  on  floppy  disks  and  subsequent  analysis. 

Two  C02  TEA  lasers  were  used  in  the  experiments:  a 
low  repetition  rate  laser  (Lumonics  model  K103),  which  can 
deliver  >  5  J/pulse  at  <0.25  Hz,  and  a  high  repetition  rate 
laser  (Tachisto,  model  555G),  which  can  deliver  <  I  J/pulse 
at  up  to  22  Hz.  For  most  experiments,  the  lasers  were  operat¬ 
ed  on  the  9.6  fi m  R  (22)  line  (1079.85  cm-  ’)  and,  in  many 
experiments,  a  three-power  Gallilean  telescope  (Janos  Tech¬ 
nology)  was  used  to  concentrate  the  beam.  Laser  pulse  ener¬ 
gy  was  monitored  continuously  with  a  volume-absorbing  ca¬ 
lorimeter  (Sctentech,  Inc.)  connected  to  a  strip-chart  record¬ 
er  Under  our  operating  conditions,  the  thyratron-triggered 
Tachisto  laser  showed  pulse-to-pulse  energy  variations  of 
about  ±5%,  while  the  older,  spark-gap-triggered  Lumon¬ 
ics  showed  variations  of  more  than  ±  20%. 

Laser  fluence  was  determined  by  the  measured  laser 
pulse  energy,  divided  by  the  laser  beam  cross-sectional  area. 
When  the  telescope  was  not  used,  the  laser  beam  cross  sec¬ 
tion  was  defined  by  an  iris  diaphram  positioned  before  the 
entrance  window.  When  the  beam  was  concentrated  by  the 
telescope,  the  area  was  determined  by  scanning  a  pinhole  in  a 
two-dimensional  survey  across  the  beam,  and  applying  the 
1/e  criterion.  The  areas  so  defined  were  very  similar  to  the 
patterns  produced  on  thermal  printer  paper,  which  was  used 
for  quick  checks  of  beam  area.  Additional  confidence  in  the 
printer-paper  measurements  was  gained  by  measuring  the 
beam  area  before,  and  after  the  telescope:  within  ±  10%, 
the  ratio  of  the  beam  areas  was  1/9.  Laser  fluence  inside  the 
cell  was  assumed  to  equal  that  measured  at  the  cell  entrance, 
because  reflection  losses  from  the  entrance  window  are  near¬ 
ly  compensated  by  back  reflections  from  the  exit  window 
and  the  small  residual  correction  is  much  less  than  the  esti¬ 
mated  +  30%  uncertainty  in  the  absolute  fluence. 

The  laser  fluence  was  varied  in  crude  steps  by  attenu¬ 
ation  with  polyethylene  films;  fine  adjustments  to  the  fluence 
were  made  by  varying  the  high  voltage  settings  on  the  laser. 
At  high  repetition  rates,  the  polyethylene  films  were  cooled 
with  a  stream  of  compressed  air  to  prevent  their  deformation 
and  melting. 

For  both  lasers,  the  far-field  pattern  had  hot  and  cold 
zones  that  differed  from  the  average  fluence  by  20%  to  30%. 
These  patterns  depended  on  laser  alignment  and  optical  path 
from  the  laser  to  the  cell,  but  they  were  somewhat  reproduc¬ 
ible,  and  their  presence  did  not  seem  to  affect  the  results.  The 
overall  intensity  profile  of  the  beam  from  each  laser  was  ap¬ 
proximately  trapezoidal,  similar  to  the  profile  measured  in 
earlier  work." 

IRMPD  yields  were  determined  using  two  variations  of 
the  VLP<t>  technique:  single-shot  accumulation,  or  fast-repe¬ 
tition  quasi-steady-state  depletion.  The  first  method  has 
been  described  previously  The  quasi-steady-state  depletion 
method  used  the  fast  repetition  rate  capability  of  the  Ta¬ 
chisto  laser  to  repeat  laser  shots  many  times  during  the  gas 
residence  time.  Since  each  parcel  of  gas  was  irradiated  many 
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times  prior  to  exiting  the  cell,  larger  depletions  of  reactant 
gas  were  generated,  allowing  measurement  of  smaller  yields- 
per-shot  that  could  be  attained  using  the  single-shot  accu¬ 
mulation  technique.  Tests  showed  that  1,1-difluoroethylene 
(the  only  ethylene  reaction  product  commercially  available: 
PCR,  Inc.)  does  not  absorb  the  laser  photons.  The  reported 
vibrational  frequencies  of  1,1-difluoroethylene  and  of  err¬ 
and  irons-  1,2-difiuoroethylene  (the  other  products)  indicate 
that  they  will  not  absorb  at  the  laser  wavelength.'3  (Yields 
were  also  measured  optoacoustically,  as  described  below.) 

In  the  absorbed  laser  energy  experiments,  the  cell  was  a 
Pyrex  cylinder  (4  cm  diam  X  30  cm  long)  equipped  with  inlet 
and  outlet  stopcocks  to  facilitate  gas  flow,  and  O-ring  con¬ 
nectors  to  accommodate  end  and  side  windows  and  the  elec¬ 
trical  feedthrough  for  the  electret  microphone/preamplifier 
iKnowles,  BT-1759).  The  output  from  the  microphone/ 
preamplifier  was  further  amplified  (Tektronix  model 
AM502)  and  captured  with  a  transient  recorder  (Biomation 
model  805)  interfaced  to  the  Nicolet  signal  averager.  By  posi¬ 
tioning  the  microphone  near  the  laser  beam,  the  initial 
acoustic  wave  generated  by  the  absorbed  energy  strikes  the 
microphone  before  the  waves  reflected  from  the  cell  walls 
can  cause  interference3;  thus,  only  the  first  maximum  in  the 
complex  signal  was  deemed  significant.  Experiments  were 
repeated  (16-256  shots!  until  a  statistical  precision  of 
±  2%-i%  was  achieved. 

The  optoacoustic  method  was  also  adapted  for  measure¬ 
ment  of  yield  vs  fluence  data  at  pressures  too  high  for  the 
VI  P<t>  technique.  The  procedure  was  to  introduce  a  known 
TFE/argon  mixture  into  the  optoacoustic  cell  and  measure 
the  OA  signal  at  a  low  fluence.  The  gas  mixture  was  then 
irradiated  for  a  known  number  of  laser  shots  at  the  higher 
fluence  to  be  investigated;  the  static  gas  mixture  was  deplet¬ 
ed  by  this  irradiation.  Following  irradiation,  the  OA  signal 
was  again  measured  at  low  fluence  to  determine  the  amount 
of  depletion,  from  which  the  IRMPD  yield-per-shot  was  de¬ 
termined,  the  reaction  products  (at  room  temperature!  do 
not  interfere  by  absorption  of  laser  light.  Small  air  leaks  li¬ 
mited  the  durations  of  data  acquisition  runs  and  the  lowest 
pressures  attainable,  but  useful  data  were  obtained,  as  de¬ 
scribed  in  the  Results  section. 

1 , 1 ,2-tnfluoroethane  was  obtained  commercially  I  PCR. 
Inc.  i  and  was  used  after  degassing  in  a  grease-free  high  vacu¬ 
um  line.  Gas  purity  was  checked  with  a  gas  chromatograph/ 
mass  spectrometer  (Varian,  Inc.  I  and  the  major  impurity  was 
found  to  be  1.1,1-tnfluoroethane  labout  0.  \%). 

RESULTS 

IRMPO  yields  v*  fluence 

The  two  versions  of  the  VLP't*  technique  were  used  to 
obtain  most  of  the  data  on  IRMPD  yields,  but  these  methods 
were  complemented  by  the  optoacoustic  technique,  which  is 
suitable  for  measuring  the  effect  of  pressure.  High  fluence 
results  were  obtained  using  the  Lumonics  laser  and  lower 
fluence  measurements  were  made  with  the  high  repetition 
rate  Tachisto  laser.  At  high  fluences,  the  pnncipal  source  of 
error  was  the  fluctuation  of  laser  power;  long-term  drift  and 
instability  of  the  mass  spectrometer  were  the  limiting  factors 
at  low  fluence.  The  minimum  uncertainty  achieved  with  the 


single-shot  accumulation  method  corresponds  to  about 
±  2%-3%,  as  determined  in  blank  runs.  Yields  of  that  or¬ 
der,  or  lower  were  not  reliably  measured  using  the  single¬ 
shot  accumulation  method  in  our  experiments.  However,  by 
using  the  high  repetition  rate  laser,  each  parcel  of  gas  enter¬ 
ing  the  VLP<I>  cell  is  irradiated  10-12  times  dunng  its  resi¬ 
dence  time,  improving  the  precision  of  the  yield  measure¬ 
ments  by  a  factor  of  5-10. 

The  results  obtained  at  1079.85  cm" 1  using  both  meth¬ 
ods  and  both  lasers  show  very  good  consistency,  even  when 
the  laser  beam  mode  structure  is  deliberately  varied.  As  a 
test,  the  area  of  irradiation  was  changed  from  2.7  to  0.78 
cm;,  using  the  Lumonics  laser,  and  the  results  were  indistin¬ 
guishable  from  those  obtained  using  a  0.42  cm;  beam  from 
the  Tachisto.  Because  of  the  insensitivity  of  the  results  to 
beam  spatial  modes,  deconvolution13  of  the  spatial  profile 
was  not  attempted,  although  an  iterative  approach  based  on 
the  calculated  results  would  be  feasible  (and  time  consum¬ 
ing). 

Tests  were  carried  out  to  eliminate  other  possible 
sources  of  experimental  artifacts.  For  example,  the  residence 
time  in  the  cell  was  held  constant  at  0.88  s  and  the  cell  pres¬ 
sure  was  varied  from  0.8  to  20  mTorr  without  significantly 
affecting  the  results  obtained;  at  higher  pressures,  the  yields 
decreased  slightly.  In  another  test,  the  gas  flow  rate  was  held 
constant  and  the  residence  time  was  varied  from  0. 10  to  0.56 
s,  without  affecting  the  results. 

All  of  these  data  are  presented  in  Fig.  1,  where  the  data 
sets  show  excellent  consistency.  The  data  are  plotted  on 
probability  graph  paper  as  a  log-normal  cumulative  distribu¬ 
tion.  because  it  has  been  shown  that  for  many  IRMPD  theo¬ 
retical  models  and  some  experimental  data,  such  a  presenta¬ 
tion  gives  straight  lines.14 

The  optoacoustic  measurements  of  yield  vs  fluence  were 
earned  out  at  pressures  ranging  from  1.5  mTorr  TFE  diluted 
in  33  mTorr  argon,  up  to  pressures  as  high  as  10  Torr,  as 
presented  in  Fig.  2.  The  most  extensive  optoacoustic  data  set 
corresponds  to  a  total  pressure  of  about  35  mTorr.  These 
data  show  the  same  fluence  dependence  as  the  collision-free 
data,  except  that  the  absolute  ytelds  are  lower,  due  to  the 
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FIG.  I  Decomposition  yields  vs  fluence  at  1071  35  cm'1  Key  O — Lu- 
monies  laser;  <} — Tachisto  laser;  • — Tachisto  laser,  very  large  number  of 
laser  shots. 
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FIG  2  Decomposition  yield  dependence  on  pressure,  determined  by  the 
optoacoustic  method.  The  curves  are  schematic.  Upper  curve:  2.7  J  cm 
absolute  uncertainties  are  ±  2%.  Lower  curve:  0. 9  J  cm  ~ absolute  uncer¬ 
tainties  are  ±  0.2%. 

higher  pressures  used  and,  perhaps,  to  small  systematic  er¬ 
rors  in  the  fluence  measurements. 

The  two  series  of  optoacoustic  measurements  presented 
in  Fig.  2  show  the  different  effects  of  pressure  on  IRMPD 
yields  at  different  laser  fluences.  At  high  fluence,  the  yield  is 
not  greatly  affected,  except  at  relatively  high  argon  pres¬ 
sures.  At  low  fluences,  however,  the  yields  are  strongly  af¬ 
fected  by  argon  pressure  over  the  entire  range  of  pressure 
investigated. 

Optoacoustic  measurements  of  absorbed  laser  power 

Calibration 

The  optoacoustic  signals  depend  on  gas  composition, 
pressure,  collisional  energy  transfer  rates,  laser  energy,  beam 
diameter,  and  spatial  position  of  the  microphone  [see  Ref. 
5(d)  for  a  fuller  discussion],  and  careful  calibration  was  nec¬ 
essary.  Absolute  calibration  factors  were  obtained  by  using 
the  absorption  cross  section  as  measured  over  a  path  length 


of  30.5  cm  and  using  the  Tachisto  laser  and  the  Scientech 
power  meter;  Beer’s  Law  fully  described  the  experimental 
data.  The  absorption  cross  section  measured  in  this  wav 
[(3.46  ±  0.35)  X  10~*  cm2]  is  about  a  factor  of  2  larger  than 
that  measured  with  a  Digilab  FTIR  with  0. 1  cm  ~ 1  resolu¬ 
tion  [(2.0  ±  0.2)  X  10" 19  cm2].  We  concluded  that  the  laser 
measurement  is  more  pertinent  to  our  experiments,  consid¬ 
ering  the  complex  spectrum  (Fig.  3)  and  possible  coupling 
mechanisms  that  may  be  operative  at  high  laser  power. 

At  sufficiently  low  fluence,  most  molecules  absorb  no 
more  than  one  photon  and  decomposition  is  not  important. 
Extensive  calibrations  were  carried  out  at  low  fluences  by 
measuring  the  microphone  signal  for  a  fixed  set  of  “hard¬ 
ware  parameters”  (e.g.,  fluence,  laser  beam  diameter,  and 
location  of  the  microphone)  and  a  fixed  partial  pressure  of 
TFE  diluted  in  various  pressures  of  argon.  The  parameters 
were  varied  systematically  to  determine  the  saturation  limit 
of  the  microphone  and  the  effects  of  the  parameters.  For 
fixed  laser  beam  spatial  parameters  and  for  TFE  highly  di¬ 
luted  in  argon,  the  microphone  signal  was  found  to  depend 
primarily  on  absorbed  laser  energy  and  slightly  on  argon 
pressure.  The  instrument  response  was  directly  proportional 
to  absorbed  laser  power  as  long  as  the  microphone  was  not 
overloaded.  The  calibration  factors  were  found  to  depend 
linearly  on  log(pressure)  for  a  wide  range  of  experimental 
parameters,  as  shown  in  Fig.  4,  where  the  calibration  factor 
is  given  by5,dl 

F  =  SoaC„/E4.  (2) 

Here,  5 OA  is  the  microphone  signal  (arbitrary  voltage  unitsi, 
Cv,  is  the  total  heat  capacity  of  the  gas  mixture  at  constant 
volume,  and  E4  is  the  absorbed  energy  per  unit  length. 
Based  on  Beer’s  law  at  the  small  absorbance  limit  and  the 
measured  cross  section  o,  E,  is  given  by 

E'=Elo[  TFE],  1 3 1 

where  EL  is  the  laser  pulse  energy  and  [TFE]  is  the  TFE 
concentration  in  convenient  units.  For  interpretation  of  the 
absorbed  laser  energy  measurements,  a  least-squares  fit  to 


FIG  3  Infrared  spectrum  of  I.I.2tnfluoroethane  near  1 100  cm  '  FTIR  spectrum  obtained  with  =0  1  cm  '  resolution,  for  vibrational  assignment,  see 
Ref  8 
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FIG  4  Optoacoustic  signal  calibration  curves.  Pressures  of  TFE  /5 — 2.4- 
24  mTorr;  ▼ — 50  mTorr  C — 125  mTorr;  solid  line:  least-squares  lit  to 
low-pressure  points. 


FIG.  5  Average  number  of  photons  absorbed  per  molecule  vs  fluence.  Ar¬ 
gon  pressures:  7—0  5  Torr,  0 — 10  Torn  /5  — 20  0  Torr,  0 — 10  0  Torr; 
dashed  line:  Master  Equation  calculation. 


the  experimental  calibration  data  was  used,  as  shown  in  the 
figure. 

Absorbed  laser  energy  measurements 

Using  the  experimental  calibration  curve,  we  can  deter¬ 
mine  the  absorbed  energy,  based  on  the  observed  micro¬ 
phone  signal  and  the  gas  composition.  To  place  the  absorbed 
energy  on  the  molecular  scale,  it  is  conveniently  expressed  as 
the  average  number  of  absorbed  photons  per  molecule,  given 
by 

(n)  =  SoaC*/|/'/iv0[TFE]),  |4) 

where  hvn  is  the  laser  photon  energy  and  the  other  factors 
have  been  identified  above.  Using  this  expression  places  the 
absorbed  energy  on  a  reduced  scale,  enabling  direct  com¬ 
parisons  of  absorbed  energies  obtained  under  widely  differ¬ 
ent  conditions. 

Expenmen ts  were  carried  out  to  determine  the  effects  of 
added  argon  and  laser  fluence  on  <n).  Usually,  better  results 
were  obtained  in  static  expcnments  than  in  flowing  gas  sys¬ 
tems,  which  suffered  from  fluctuations  in  total  pressure  and 
in  gas  composition.  Care  was  taken  to  avoid  depletion  of  the 
TFE  in  experiments  with  high  laser  fluence.  The  results  are 
presented  in  Fig.  5,  and  they  show  no  significant  dependence 
on  pressure  for  fluences  up  to  about  2  J  cm  " 2.  The  results  are 
highly  consistent,  with  a  relative  uncertainty  of  about  ±  5 % 
and  an  absolute  uncertainty  estimated  at  about  ±  30%, 
mostly  due  to  the  uncertainty  in  the  fluence  and  energy  mea¬ 
surements  that  went  into  calibration  of  the  microphone. 

At  the  highest  fluences,  decomposition  of  the  TFE  dur¬ 
ing  the  laser  pulse  may  become  significant,  affecting  the  sim¬ 
ple  interpretation  of  the  results.  Although  the  reactions  are 
slightly  endothermic  and  can  affect  the  optoacoustic  signal 
to  a  small  extent,  the  major  uncertainty  is  due  to  the  un¬ 
known  extent  of  absorption  of  the  laser  photons  by  vibra¬ 
tionally  hot  reaction  products.  In  the  infrared  fluorescence 
measurements  on  this  system/  direct  evidence  was  obtained 
for  secondary  photolysis  of  the  products. 


MASTER  EQUATION  MODEL 

The  basic  Master  Equation  approach  to  simulating  ex¬ 
periments  such  as  these  has  been  described  in  detail,6  except 
for  several  significant  improvements.  Essentially,  the  ap¬ 
proach15  used  Monte  Carlo  techniques  to  simulate  random 
walks  of  individual  molecules  in  energy  space,  subject  to  the 
properties  of  the  molecule  and  its  interactions  with  the  laser 
photons.16  Many  such  “trajectories”  are  calculated  and  the 
results  are  averaged,  or  otherwise  stored  for  examination. 
Other  than  that  used  in  calculating  the  vibrational  densities 
of  states,  no  intrinsic  “graining”  is  imposed  on  the  physics  of 
the  system,  unlike  approaches  that  use  an  “energy-grained” 
Master  Equation.  Thus,  convergence  of  the  solutions  does 
not  depend  on  grain  size,  although  the  precision  obtainable 
depends  on  the  number  of  trajectories  calculated,  as  is  usual 
in  Monte  Carlo  simulations.16 

The  physical  processes  considered  in  the  model  include 
optical  absorption  and  stimulated  emission,  collisional  ener¬ 
gy  transfer,  infrared  fluorescence,  and  up  to  three  unimole- 
cular  reaction  channels;  several  models  for  some  of  these 
processes  have  been  incorporated  in  the  computer  code  as 
options.6 

The  Master  Equation  implementation  has  been  signifi¬ 
cantly  improved  in  several  ways  for  present  purposes.  First, 
all  densities  of  vibrational  states  are  now  calculated  from 
exact  counts  of  states,17  rather  than  from  the  Whitten-Ra- 
binovitch  approximation18  used  earlier.  For  TFE,  an  energy 
grain  of  25  cm  " 1  has  proved  to  be  adequate  for  the  calcula¬ 
tions.  The  densities  of  states  for  the  first  2500  cm  - 1  of  inter¬ 
nal  energy  are  stored  in  a  100-element  array;  in  addition,  the 
densities  for  the  first  50  000  cm  * 1  are  stored  in  a  second  100- 
element  array.  The  density  of  states  at  any  energy  up  to 
50  000  cm  ~ 1  is  selected  by  interpolation  between  elements  in 
one  of  these  two  arrays.  Incorporation  of  exact  state  counts 
has  made  little  difference  in  computed  results,  except  for  the 
distribution  of  initial  energies  selected.  For  example,  now 
the  thermal  distribution  is  much  more  accurately  produced, 
compared  to  the  somewhat  shifted  version  obtained6  using 
the  Whitten-Rabinovitch  approximation. 
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The  other  improvements  are  concerned  with  the  optical 
pumping  and  stimulated  emission  segments  of  the  model.  In 
addition  to  the  former  capability  to  consider  a  laser  intensity 
that  is  constant  for  the  entire  duration  of  a  trajectory  and  a 
laser  pulse  shape  that  decays  exponentially,  additional  pro¬ 
vision  has  been  made  for  a  laser  pulse  of  constant  duration 
and  finite  length.  This  capability  enables  investigation  of 
post-pulse  phenomena  subsequent  to  a  pulse  of  constant  in¬ 
tensity.  The  second  enhancement  of  the  optical  interaction 
model  was  to  incorporate  Quack’s  theory  of  infrared  laser 
absorption/emission  for  his  cases  B  and  C. 19 

Prior  to  incorporation  of  Quack’s  theory  in  our  model, 
we  used  a  semiempirical  expression  for  the  absorption  cross 
section  dependence  on  vibrational  energy.4-61416 

c*£)  =  <70[1  +£7Avora'.  (5) 


where  a0  is  the  absorption  cross  section  for  all  molecules  at 
the  zero-point  energy  level,  hv0  is  the  laser  photon  energy, 
and  «U  I  is  a  parameter  that  depends  on  the  wavelength  of  the 
laser.  This  expression  has  been  used  in  Master  Equation 
treatments  of  IRMPD  with  fair  success,  but  recently,  Haas 
and  co-workers  commented  on  the  inadequacies  of  a  Master 
Equation  treatment  that  employed  this  type  of  expression.2 
In  order  to  reduce  the  degree  of  empiricism  somewhat,  and 
because  we  also  find  Eq.  ( 5 )  to  be  inadequate  for  treatment  of 
the  results  obtained  in  the  present  work  (see  below),  we  im¬ 
plemented  the  Quack  formalism  for  cases  B  and  C. 

Quack  has  described  his  statistical-dynamical  theory 
extensively  in  the  literature19  and  we  will  not  attempt  to  re¬ 
peat  the  details  here.  His  approach  is  to  evaluate  the  matrix 
elements  for  laser  pumping  and  stimulated  emission  using  a 
combination  of  ab  initio  theory  and  judicious  approxima¬ 
tions.  The  resulting  expressions  give  the  optical  transition 
rates  in  terms  of  physically  measurable  quantities,  with  pro¬ 
vision  for  empirical  parametrization,  where  desirable. 

According  to  Quack,'9  the  optical  transition  rates  for 
up-pumping  from  level  M  to  M  4-  1  for  both  cases  B  and  C 
are  given  by  the  expression 

K w  -  i.,w  =  _  i  v/ 1  ■/$.*_  i ,  |6) 

where  „  is  the  frequency  spacing  of  effectively  cou¬ 
pled  states  in  the  upper  level  and  I  V%,  *  ,  M  is  the  absolute 
value  of  the  matrix  element  connecting  the  two  levels  leach 
microcanomcal  level  contains  many  individual  quantum 
states).  According  to  Quack,  8,  can  be  related  to  the  total 
density  of  vibrational  states  according  to  the  expression 

<5,  =  2irc/p\E,),  }7| 


where  c  is  the  speed  of  light,  andpi£,  I  is  the  total  density  of 
v  ibrational  states  at  the  energy  of  level  i.  The  matrix  element 
is  related  to  the  laser  intensity  /  and  the  integrated  absorp¬ 
tion  cross  section  |  j  odx  | ,  by  the  expression 
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where  p'<E)  is  the  density  of  vibrational  states  with  the 
pumped  vibrational  mode  in  the  zero  level,  1EQ  is  the  level 
width  over  which  the  matrix  elements  are  approximately 
constant,  and  hvp  is  the  energy  (expressed  in  cm  “ ')  of  the 
pumped  mode  The  factor  5  £w  i  allows  empirical  parame- 


tnzation  of  the  expression;  when  £(£„)=  i  .0,  the  param¬ 
eters  are  assumed  to  be  independent  of  internal  energy  and  to 
equal  those  measured  at  low  intensities.  For  intensities  ex¬ 
pressed  in  W  cm-2,  state  densities  in  states  per  cm' ',  a  in 
cm2,  A£e  in  cm-1,  and  dx  in  cm-',  the  proportionality 
constant  is  a  =  1.509  X  1033,  to  express  the  squared  matrix 
element  in  units  of  rad2  s-2. 

The  down-pumping  rates  differ  for  cases  B  and  C.  For 
case  B,  the  down-transition  rate  from  level  M  to  M  -  1  is 
given  by 
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where  the  factors  have  been  defined  above.  For  case  C,  the 
down-transition  rate  is  given  by 
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where  the  symbols  have  the  same  meanings  and  units  as 
above,  and  the  proportionality  constant  y  is  1 . 784  x  1 0  ~ 6  for 
the  transition  rate  expressed  in  s _ The  crossover  from  case 
C  to  case  B  at  low  energy  and/or  low  intensity  is  treated  by 
selecting  the  larger  of  the  two  calculated  transition  rates,  as 
described  by  Quack. 

For  1,1,2-tnfluoroethane,  the  molecular  parameters 
needed  for  use  of  the  Quack  expressions  are  known  and  are 
summarized  in  Tables  1  and  II.  The  expressions  can  be  em¬ 
pirically  adjusted19  by  treating  A EQ  as  a  parameter  used  in 
concert  with  5  l£  I,  as  discussed  below.  Densities  of  vibration¬ 
al  states  were  obtained  from  exact  counts  by  using  the  vibra¬ 
tional  assignment  for  the  molecule  and  assuming  harmonic 
oscillators.  For  the  density  of  states  calculations,  the  internal 
rotor  was  treated  as  a  1 17  cm  "  1  vibration. 


DISCUSSION 

The  two  experimental  techniques  employed  in  this  in¬ 
vestigation  provide  information  about  the  IRMPA  process 
over  two  complementary  internal  energy  ranges.  The  optoa- 
coustic  measurements  are  most  accurate  for  low  fluences,  or 
when  reaction  is  not  important.  On  the  other  hand,  the 
IRMPD  measurements  require  that  molecules  be  pumped 
above  the  reaction  threshold,  and  it  has  been  shown  that  the' 
technique  provides  information  mostly  about  the  reaction 
threshold  region  and  above.4 14 16 

The  absorbed  energy  measurements  show  little  depen¬ 
dence  on  pressure,  and  thus  we  conclude  that  at  low  pres¬ 
sures.  where  collisional  deactivation  during  the  pumping 
process  is  unimportant,  the  measurements  reflect  the  aver¬ 
age  energy  of  the  nascent  population.  Since  the  reaction 
threshold  energies  are  about  24  000  cm-1,  the  optoacoustic 
measurements  can  be  interpreted  in  a  straightforward  way 
up  to  (n)  is  22.  Above  that  value,  the  extent  of  decomposi¬ 
tion  will  affect  the  results  in  three  ways:  1 1 1  depletion  of  the 
parent  molecules  during  the  laser  pulse  will  mean  fewer  pho¬ 
tons  absorbed,  i2l  the  reaction  thermochemistry  can  affect 
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TABLE  I.  Molecular  properties  of  1.1.2-tnfluoroethane. 


Property 

Value 

Reference 

Vibrational  frequency  assignment 

a 

Heat  capacity  |C.) 

16.29  *■  0.033  ( T  — 

3001  Gibbs  mol ' 1 

b 

Lennard-Jones  a 

5  28  A 

C 

Lennard-Jones  e/k 

325  K 

hvf  i  pumped  mode:  v,,| 

10?6  cm  ~ ' 

ja  dx 

i5  2:  1)  x  10"  '*  cm 

d 

L'nunolecuiar  reaction  parameters 

1  Whitten-Rabuiovi  tch ) 

e 

Reaction:  A 

B 

c 

<^v>TTT 

941.62 

917  37 

906.77 

£.m 

23  714 

23  714 

24413 

TST  Zero  point  energy 

9  751.0 

9  538.5 

9  64  1  0 

Whitten-Rabmovitch  0 

1.353 

1  377 

1  379 

<'/r  ft) 

0  65 

2.50 

2  02 

*V  F.  Kalasinsky,  H.  V.  Aryans,  and  T.  S.  Little,  J.  Phys.  Chem.  86,  1351 1 1982). 

'Calculated  from  the  vibrational  assignment. 

'  Estimated  using  Lyderson's  method  as  described  in  R.  C  Reid  and  T  1C  Sherwood,  The  Properties  of  Gases 
and  Liquids,  2nd  ed.  iMcGraw-Hill,  New  York,  1966). 

1  Measured  in  this  work;  see  FTIR  spectrum  in  Fig.  3. 

•  Whitten-Rabmovitch  parameters  for  the  transition  states  were  denved  from  the  detailed  transition  state 
frequencies  given  by  B.  E.  Holmes.  D.  W.  Setser.  and  G.  O.  Pritchard,  Int  J  Chem.  Kinet.  8.  215  1 19761; 
RRK.M  k \E  i's  were  calculated  using  the  TFE  state  densities  and  Whitten-Rabmovitch  approximation  for  the 
sums  of  states. 


the  heat  release,  and  (3)  the  [hot]  reaction  products  produced 
during  the  pulse,  or  in  previous  pulses,  may  absorb  the  laser 
light.  Inspection  of  Fig.  5  shows  that  with  10  Torr  argon 
added,  (n)> 32,  much  higher  than  the  reaction  threshold 
energy.  At  this  high  pressure,  decomposition  is  inhibited  by 
collisional  deactivation  of  excited  molecules  and  the  reac¬ 
tant  is  not  depleted.  Indeed,  the  collisional  deactivation  rate 
is  fast  enough  so  that  the  excited  molecules  are  never 
pumped  very  high  on  the  energy  ladder. 

The  absence  of  a  pressure  dependence  at  relatively  high 
pressures,  where  collisional  deactivation  during  the  laser 
pulse  is  important,  indicates  that  the  net  optical  absorption 
rate  of  the  excited  molecules  does  not  depend  strongly  on 
internal  energy.  Inspection  of  the  10  Torr  added  argon  data 
set  in  Fig.  5  shows  that  it  is  virtually  indistinguishable  from 
data  sets  obtained  at  much  lower  pressures.  More  quantita¬ 
tive  conclusions  about  the  optical  pumping  process  require 
the  numerical  calculations  described  below. 

The  IRMPD  measurements  probe  molecules  pumped 
above  the  reaction  threshold.  When  the  yields  are  low,  the 
decomposing  molecules  represent  only  the  high  energy 
tail"  of  the  population  distribution,  which  is  highly  sensi¬ 
tive  to  several  effects,  including  the  optical  pumping  process 
and  collisional  energy  transfer.  At  sufficiently  low  pressures, 
homogeneous  collisional  effects  are  unimportant,  and  the 
IRMPD  yields  reflect  the  fraction  of  molecules  pumped 
above  the  reaction  threshold  that  react  before  deactivation  at 


TABLE  II  Quack  theory  lined  parameters 


Parameter 

Value 

Pulse  width 

125  £  15  ns 

XEq 

4  £  1  cm  ' 

w 

23  700  ;  100  cm  '' 

m 

20  £  5 

the  walls  of  the  cell.  Since  drift  times  to  the  cell  walls  are  of 
the  order  of  3  x  10"4  s,  the  decomposition  yield  at  low  pres¬ 
sure  equals  the  fraction  of  molecules  that  attain  energies 
high  enough  to  have  specific  unimolecular  rate  constants 
&l£)>3x  10J  s_1  (i.e.,  £>26  000 cm" '). 

Information  about  the  optical  pumping  process  can  be 
obtained  from  the  dependence  of  the  collision-free  IRMPD 
yield  on  fluence.  As  a  class,  systems  with  optical  pumping 
properties  that  can  be  described  by  the  simple  power-law 
equation  1 51  have  decomposition  yields  that  can  be  plotted 
linearly  vs  logifluencel  on  probability  graph  paper,  as  has 
been  discussed  elsewhere. 14  16  The  VLP4>  experimental  data 
are  presented  on  such  a  log-normal  plot  in  Fig.  1  and  the 
resulting  line  is  clearly  curved,  indicating  that  the  optical 
pumping  is  not  well  described  by  Eq.  15).  The  heuristic  ap¬ 
proach  taken  in  our  earlier  work  in  order  to  account  for  such 
curvature  was  to  assume  that  the  excited  molecules  were 
partitioned  into  two  or  more  population  subsets  with  differ¬ 
ent  optical  pumping  parameters.'0  This  approach  could 
have  been  taken  in  the  present  work,  but  Quack's  more  so¬ 
phisticated  theory"*  accounts  satisfactorily  for  the  expen- 
mental  observations  without  recourse  to  the  multiple-popu¬ 
lation  model. 

A  more  complete  discussion  of  collisional  energy  trans¬ 
fer  tn  this  system  will  appear  elsewhere,:i  but  IRMPD  yield 
data  obtained  with  the  optoacoustic  technique  are  presented 
in  Fig.  2.  These  results  show  that,  as  expected,  the  low 
fluence  data  are  far  more  susceptible  to  collisional  effects 
than  data  obtained  at  higher  fluences.  The  high  energy  tail  of 
the  excited  molecule  distnbution  is  sensitively  affected  by 
collisions,  even  at  pressures  as  low  as  35  mTorr.  iNote  that 
our  expenmental  tests  showed  that  the  VLP4>  data  obtained 
at  20  mTorr  were  not  significantly  affected  by  total  pres¬ 
sure!  The  decreased  yield  is  due  almost  entirely  to  "post¬ 
pulse"  collisional  deactivation  of  excited  molecules  that 
have  energies  only  a  little  above  the  reaction  threshold.  In 
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addition,  mass  transport  to  the  cell  wall  is  affected  by  the 
higher  pressures,  and  the  detailed  differences  in  the  collision 
efficiency  of  walls  vs  gas  phase  must  be  taken  into  account 
for  a  full  description. Such  a  description  is  feasible,  but  it  is 
beyond  the  scope  of  the  present  work. 

To  obtain  more  quantitative  information  about  the 
pumping  process,  the  Master  Equation  approach  described 
above  was  utilized  and  the  effects  of  many  parameters  were 
investigated.  Insofar  as  possible,  experimentally  measured 
parameters  were  used  in  order  to  limit  the  latitude  of  the 
calculations.  Many  combinations  of  assumptions  and  input 
parameters  were  investigated  to  determine  whether  they 
were  consistent  with  the  experimental  data,  but  it  was  found 
that  the  two  sets  of  experimental  data  combined  with  knowl¬ 
edge  of  TFE  molecular  properties  constrained  the  latitude  of 
the  calculations  greatly,  and  a  satisfactory  simulation  of  the 
experiments  was  difficult  to  achieve. 

The  unsatisfactory  attempts  at  simulation  of  the  experi¬ 
ments  included  the  following  combinations  of  assumptions: 

1 1 1  Exponential  or  square  laser  pulse  used  with  the  pow¬ 
er-law  absorption  cross  section  in  Eq.  (5):  This  combination 
could  be  made  to  fit  the  absorbed  energy  measurements  very 
well  by  using  the  measured  cross  section  and  adjusting  the 
parameter  n\A  ),  but  the  predicted  yield  vs  fluence  results 
were  too  steep  on  a  log-normal  plot,  as  shown  in  our  prelimi¬ 
nary  report.'3  Also,  the  simulated  log-normal  plot  was  a 
straight  line,  unlike  the  experimental  data. 

'2|  Exponential  laser  pulse  and  Quack's  theory  with 
5  :£  I  =  1.0:  Here,  the  absorbed  energy  measurements  could 
be  well  simulated  by  varying  the  parameter  A £c,  but  the 
predicted  reaction  yields  were  much  too  high.  Furthermore, 
the  pulse  decay  time  significantly  affected  the  predicted  val¬ 
ue  of  (n)  at  a  given  fluence.  This  is  because  the  transition 
between  cases  B  and  C  depends  on  laser  intensity,  and  the 
down-pumping  rate  in  case  C  is  much  greater  than  in  case  B. 
Molecules  pumped  up,  early  in  the  pulse  under  case  B  condi¬ 
tions.  could  be  rapidly  down  pumped  later  in  the  pulse,  when 
the  intensity  is  lower  and  case  C  is  operative.  This  effect  may 
be  observable  in  carefully  designed  experiments;  its  presence 
may  indicate  that  an  accurate  simulation  of  the  laser  pulse  is 
necessary  for  more  detailed  calculations. 

1 3 1  Square  laser  pulse  and  Quack's  theory  with 
c  £  i  =  1.0:  To  simplify  the  model,  a  square  pulse  was  as¬ 
sumed.  This  simplified  version  gave  good  agreement  with 
< n )  for  suitable  choices  of  A£c  and  the  pulse  width,  but  the 
predicted  reaction  yields  were  much  too  high,  indicating 
that  the  average  molecules  couid  be  well  modeled,  but  the 
high-energy  tail  of  the  distribution  Ithat  resulted  in  decom¬ 
position!  was  predicted  to  be  too  highly  populated. 

Since  the  Quack  theory  with  £(£)  =  10  land  appropri¬ 
ate  choices  of  XEq  and  laser  pulse  widthi  fitted  the  experi¬ 
mental  data  on  absorbed  laser  energy  very  well,  our  ap¬ 
proach  was  to  choose  a  functional  form  that  gave  51  £1  =  10 
at  low  internal  energies,  but  allowed  simple  variations  at 
higher  energies.  A  first  attempt  at  a  functional  form  for  1 1£  1 
incorporated  a  Lorentzian  with  “width”  and  “center”  that 
could  vary  linearly  as  a  function  of  internal  energy,  as  sug¬ 
gested  by  Quack. 'g'dl  This  approach  was  not  successful,  al¬ 
though  nonlinear  dependencies  on  internal  energy  could 


probably  be  made  to  work.  Inspection  of  the  infrared  spec¬ 
trum  of  TFE  (Fig.  3)  shows  that  at  least  two  fundamental 
transitions  are  important,  and  a  simple  Lorentzian  is  not 
adequate. 

The  second  attempt  at  choosing  a  functional  form  for 
5  l£)  was  successful.  The  assumed  form  incorporates  only 
two  empirical  parameters  and  it  gives  £  1  £ )  =  1.0  at  low  ener¬ 
gies: 

$  (£1  =  exp|  -  \E/W r|.  (12| 

Although  this  function  is  arbitrary  and  not  unique,  it  pro¬ 
vides  a  smooth  transition  between  the  optical  properties  of 
the  molecule  at  low  energies,  where  they  can  be  measured 
with  a  spectrometer,  and  at  high  energies,  where  the  anhar- 
monicities  may  become  very  significant. 

By  a  process  of  trial  and  error,  a  combination  of  param¬ 
eters  (Table  II)  was  found  that  gives  model  simulations  in 
good  agreement  with  the  experimental  data  sets,  as  shown  in 
Figs.  1  and  5.  The  calculated  yield  vs  fluence  results  are  in 
excellent  agreement  with  the  experimental  data  over  virtual¬ 
ly  the  entire  fluence  range  investigated.  The  minor  discrep¬ 
ancy  at  very  low  fluence  may  be  due  to  experimental  artifacts 
and/or  to  inadequacies  of  the  model  for  the  extremely  small 
populations  at  the  highest  part  of  the  distribution.  The  calcu- 


cZ  0 


FIG  6  Population  distributions  Each  panel  corresponds  to  a  different  13- 
'er  fluence  For  reference,  the  300  K  thermal  distnbution  and  the  function 
c  £  'jreshowninthetoppanelithepopulaiiondistnbutionfor  I  5J  cm  is 
essentially  zero  below  5000  cirr'i 
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TABLE  I  Molecular  properties  of  l.!.2-trifluoroethane. 


Property 

Value 

Reference 

Vibrational  frequency  assignment 

a 

Heel  capacity  (C, ) 

16.29  +  0.033  IF- 

300)  Gibbs  mol-' 

b 

Lennard-Jones  o 

5.28  A 

c 

Lennard-Jones  e/k 

325  K 

hvr  ipumped  mode:  vM| 

1076  cm-1 

So  dx 

|5±  llXl0-"cm 

d 

Unimolecular  reaction  parameters 

( Whitten-Rabinovitch  | 

e 

Reaction:  A 

B 

c 

(hv)rvT 

941.62 

917  37 

906.77 

23  714 

23  714 

24  413 

TST  Zero  point  energy 

9  751  0 

9  538.5 

9641  0 

Whitten-Rabinovitch  0 

1.353 

1  377 

I  379 

<V*xt  /*//> 

0.65 

2  50 

2.02 

*  V.  F.  Kalasinsky,  H.  V.  Anjana,  and  T.  S.  Little.  !.  Phys.  Chem.  86,  1351  ( 1 982). 

"Calculated  from  the  vibrational  assignment 

■  Estimated  using  Lyderson's  method  as  descnbed  in  R.  C.  Reid  and  T.  K.  Sherwood,  The  Properties  of  Oases 
and  Liquids,  2nd  ed.  (McGraw-Hill,  New  York.  19661. 

“  Measured  in  this  work:  see  FTIR  spectrum  in  Fig.  3. 

•  WHitten-Rabinovitch  parameters  for  the  transition  states  were  denved  from  the  detailed  transition  state 
frequencies  given  by  B.  E.  Holmes,  D.  W.  Setser.  and  G.  O.  Pritchard.  Int.  J.  Chem.  Kinet.  8.  215  (1976); 
RRKM  k \E  |'s  were  calculated  usmg  the  TFE  state  densities  and  Whitten-Rabtnovitch  approximation  for  the 
sums  of  states. 


the  heat  release,  and  (3)  the  [hot]  reaction  products  produced 
during  the  pulse,  or  in  previous  pulses,  may  absorb  the  laser 
light.  Inspection  of  Fig.  5  shows  that  with  10  Torr  argon 
added.  (n>> 32,  much  higher  than  the  reaction  threshold 
energy.  At  this  high  pressure,  decomposition  is  inhibited  by 
coiltsional  deactivation  of  excited  molecules  and  the  reac* 
tant  is  not  depleted.  Indeed,  the  collisional  deactivation  rate 
is  fast  enough  so  that  the  excited  molecules  are  never 
pumped  very  high  on  the  energy  ladder. 

The  absence  of  a  pressure  dependence  at  relatively  high 
pressures,  where  collisional  deactivation  during  the  laser 
pulse  is  important,  indicates  that  the  net  optical  absorption 
rate  of  the  excited  molecules  does  not  depend  strongly  on 
internal  energy.  Inspection  of  the  10  Torr  added  argon  data 
set  in  Fig.  5  shows  that  it  is  virtually  indistinguishable  from 
data  sets  obtained  at  much  lower  pressures.  More  quantita¬ 
tive  conclusions  about  the  optical  pumping  process  require 
the  numerical  calculations  descnbed  below. 

The  IRMPD  measurements  probe  molecules  pumped 
above  the  reaction  threshold.  When  the  yields  are  low,  the 
decomposing  molecules  represent  only  the  high  energy 
“tail”  of  the  population  distribution,  which  is  highly  sensi¬ 
tive  to  several  effects,  including  the  optical  pumping  process 
and  collisional  energy  transfer.  At  sufficiently  low  pressures, 
homogeneous  collisional  effects  are  unimportant,  and  the 
IRMPD  yields  reflect  the  fraction  of  molecules  pumped 
above  the  reaction  threshold  that  react  before  deactivation  at 


TABLE  II.  Quack  theory  fitted  parameters. 


Parameter 

Value 

Pulse  width 

125  £  15  ns 

4  +  1  cm  ‘ 1 

w 

23  700  ±  100  cm  '' 

m 

20  ±  5 

the  walls  of  the  cell.  Since  drift  times  to  the  cell  walls  are  of 
the  order  of  3  X  10-4  s,  the  decomposition  yield  at  low  pres¬ 
sure  equals  the  fraction  of  molecules  that  attain  energies 
high  enough  to  have  specific  unimolecular  rate  constants 
A(£)>3x  103  s~‘  (i.e.,  £>26  000  cm  ~‘). 

Information  about  the  optical  pumping  process  can  be 
obtained  from  the  dependence  of  the  collision-free  IRMPD 
yield  on  fluence.  As  a  class,  systems  with  optical  pumping 
properties  that  can  be  descnbed  by  the  simple  power-law 
equation  ( 5)  have  decomposition  yields  that  can  be  plotted 
linearly  vs  log(fluence|  on  probability  graph  paper,  as  has 
been  discussed  elsewhere. 14 16  The  VLP<t>  expenmen tal  data 
are  presented  on  such  a  log-normal  plot  in  Fig.  1  and  the 
resulting  line  is  clearly  curved,  indicating  that  the  optical 
pumping  is  not  well  described  by  Eq.  1 5).  The  heuristic  ap¬ 
proach  taken  in  our  earlier  work  in  order  to  account  for  such 
curvature  was  to  assume  that  the  excited  molecules  were 
partitioned  into  two  or  more  population  subsets  with  differ¬ 
ent  optical  pumping  parameters.'0  This  approach  could 
have  been  taken  in  the  present  work,  but  Quack's  more  so¬ 
phisticated  theory19  accounts  satisfactorily  for  the  experi¬ 
mental  observations  without  recourse  to  the  multiple-popu¬ 
lation  model. 

A  more  complete  discussion  of  collisional  energy  trans¬ 
fer  in  this  system  will  appear  elsewhere,21  but  IRMPD  yield 
data  obtained  with  the  optoacoustic  technique  are  presented 
in  Fig.  2.  These  results  show  that,  as  expected,  the  low 
fluence  data  are  far  more  susceptible  to  collisional  effects 
than  data  obtained  at  higher  fluences.  The  high  energy  tail  of 
the  excited  molecule  distribution  is  sensitively  affected  by 
collisions,  even  at  pressures  as  low  as  35  mTorr.  iNote  that 
our  experimental  tests  showed  that  the  VLP<P  data  obtained 
at  20  mTorr  were  not  significantly  affected  by  total  pres¬ 
sure.!  The  decreased  yield  is  due  almost  entirely  to  "post¬ 
pulse"  collisional  deactivation  of  excited  molecules  that 
have  energies  only  a  little  above  the  reaction  threshold.  In 
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lated  curve  for  (n)  vs  fluence  in  Fig.  5  corresponds  to  an 
assumed  argon  pressure  of  1 .0  Torr,  and  the  saturation  effect 
due  to  the  function  § (E )  is  clearly  seen;  at  higher  pressures, 
molecules  are  deactivated  before  reaching  such  high  ener¬ 
gies.  and  the  predicted  behavior  then  resembles  that  ob¬ 
served  in  the  experiments  with  10  Torr  argon. 

The  fitted  parameters  appear  to  be  reasonable,  but  they 
are  not  unique,  because  they  depend  on  the  assumed  func¬ 
tional  form  for  $  (£  1.  Within  the  framework  of  our  assump¬ 
tions,  however,  there  is  very  little  latitude  for  variation  of  the 
parameters,  as  indicated  by  the  estimated  uncertainties. 

Comparison  of  the  parameters  in  Table  II  with  the 
known  molecular  properties  of  TFE  and  the  experimental 
conditions  shows  good  agreement.  For  example,  the  laser 
pulse  width  of  125  ±  15  ns  for  a  simulated  square  pulse  is 
roughly  similar  to  the  typical  CO,  TEA  laser  pulse,  which  is 
comprised  of  many  irratically  modulated  spikes  enclosed  by 
an  average  amplitude  envelope  consisting  of  a  100  ns  pulse 
with  a  1-2 /is  tail.  Also,  the  parameter  A EQ  is  expected  to  be 
of  the  same  order  as  the  width  of  the  infrared  absorption 
band19;  in  TFE,  the  laser  photons  are  4  cm  - 1  out  of  reso¬ 
nance  with  the  vn  Q  branch  (for  the  C,  conformer,  the  pre¬ 
dominant  form  in  the  gas  phase)  at  1076  cm-1,  which  is 
about  2. 5  cm - 1  in  width,  compared  to  XEQ  =  4  +  1  cm'1. 
Most  interesting,  however,  is  the  fitted  value  of  the  param¬ 
eter  W. 

In  Eq.  (12),  the  parameter  W  is  the  energy  at  which 
S  I  £ )  =  l/e  and  the  optical  coupling  matrix  element  is  much 
reduced  compared  to  its  value  for  transitions  at  low  internal 
energies.  There  can  be  several  reasons  for  a  reduction  in 
magnitude  of  the  matrix  elements,  but  a  good  explanation  is 
provided  if  the  transition  energies  shift  away  from  resonance 
due  to  anharmonicity.  The  value  of  W  depends  slightly  on 
that  of  the  parameter  m,  but  for  m  =  20,  the  uncertainty  in 
W  is  quite  small. 

The  fact  that  the  value  of  W  is  so  close  to  the  reaction 
thresholds  is  highly  suggestive  that  anharmonicity  intro¬ 
duced  by  the  open  reaction  channels  plays  an  important  role 
in  the  spectroscopic  properties  of  this  highly  vibrationaiiy 
excited  molecule.  Although  §  i£ )  is  not  a  unique  function,  its 
general  behavior  may  provide  a  good  semiquantitative  de¬ 
scription  of  the  optical  matrix  elements  near  the  reaction 
thresholds.  Based  upon  our  experience  in  attempting  to  sim¬ 
ulate  the  experimental  data  presented  here,  we  conclude 
that,  within  the  framework  of  the  Quack  theory,  the  matrix 
elements  must  undergo  a  drastic  change  near  the  reaction 
thresholds.  If  our  conclusion  is  correct,  this  may  be  the  first 
time  the  effect  has  been  observed  in  highly  vibrationaiiy  ex¬ 
cited  polyatomic  molecules. 

The  reactions  in  TFE  are  HF  eliminations  that  must 
pass  through  molecular  configurations  starting  from  single- 
bonded  carbon  atoms  with  sp}  hybridization,  through  inter¬ 
mediate  configurations  ii.e„  a  transition  statei,  to  doubly 
bonded  carbon  atoms  with  sp:  hybridization.  The  vM  in¬ 
frared  transition  at  1076  cm'1  is  an  asymmetric  C-F 
stretching  mode,  a  mode  that  must  be  affected  by  a  change  in 
hybridization  of  the  carbon  atom:  as  the  hybridization 
changes  from  sps  to  spl,  the  C-F  bond  will  be  slightly 


cies.  Specifically,  the  1076  cm  ~ 1  C-F  stretch  mode  frequen¬ 
cy  in  TFE  may  be  compared  to  the  C-F  stretch  modes  in  cis- 
and  rrans-l,2,-difluoroethylene  (product  molecules),  which 
have  frequencies  of  1014, 1127, 1123, and  1 159cm~ ',  far  out 
of  resonance  with  the  laser  photons.  Since  every  C-F  bond  in 
TFE  (even  those  not  directly  involved  in  the  HF  elimination) 
undergoes  a  change  in  hybridization  when  reaction  occurs,  it 
is  not  surprising  that  the  transition  matrix  elements  for  the 
C-F  stretch  modes  drastically  decrease  in  magnitude  at  en¬ 
ergies  close  to  the  reaction  thresholds. 

The  primary  objective  of  this  work  was  to  determine  the 
population  distributions  produced  in  TFE  by  IRMPA.  Al¬ 
though  not  unique,  the  model  results  probably  represent  a 
very  good  approximation  to  the  distribution  functions,  espe¬ 
cially  considering  the  narrow  latitudes  available  for  satisfac¬ 
tory  simulation  of  the  two  data  sets.  The  population  distribu¬ 
tions  predicted  for  several  laser  fluences  are  presented  in  Fig. 
6  along  with  a  plot  of£  (£ ).  Note  that  the  250cm  “ 1  “binning" 
of  the  calculated  results  was  for  bookkeeping  purposes  and 
had  no  effect  on  the  actual  calculation. 

The  calculated  distributions  predict  that  at  low  and  in¬ 
termediate  fluences,  a  significant  fraction  of  molecules  re¬ 
mains  m  the  initial  thermal  distribution,  but  that  fraction 
becomes  much  smaller  at  higher  fluences.  At  all  fluences,  the 
initial  thermal  distribution  is  translated  upwards  by  integral 
numbers  of  photon  energies,  resulting  in  a  series  of  peaks 
that  have  about  the  same  width  as  the  initial  thermal  distri¬ 
bution.  Also,  the  overall  dispersion  of  energies  is  quite  large, 
compared  to  the  mean.  At  high  fluences,  however,  the  dis¬ 
persion  becomes  much  smaller,  as  molecules  “stack  up"  at 
high  energies,  due  to  the  fall  off  of  £  (£).  On  a  longer  time 
scale,  most  molecules  will  decompose,  because  they  have  en¬ 
ergies  greater  than  the  reaction  thresholds. 

These  population  distributions  illustrate  that  IRMPA  is 
effective  in  exciting  most  molecules  in  the  laser  beam  and 
that  the  distributions  are  quite  broad.  For  effective  use  of  the 
technique  in  laboratory  studies  of  other  molecules,  efforts 
are  necessary  to  elucidate  the  distribution  functions  pro¬ 
duced.  Absorbed  energy  measurements  taken  alone  are  not 
sufficient  to  determine  the  population  distributions,  since 
any  of  the  combinations  of  assumptions  described  above 
could  simulate  the  cptoacoustic  experiments. 

The  combinations  of  experimental  measurements  that 
are  most  effective  in  constraining  the  computer  simulations 
are  those  that  measure  different  attributes  of  the  population 
distributions.  For  example,  in  the  present  study,  the  ab¬ 
sorbed  energy  measurements  gave  information  mostly  about 
the  mean  energy  ithe  first  moment  of  the  population  distri¬ 
bution),  while  the  decomposition  yield  vs  fluence  measure¬ 
ments  probed  the  high  energy  tail  of  the  distribution  func¬ 
tion.  In  Fig.  7,  <£)  is  plotted  as  a  function  of  (n)  toillustrate 
the  direct  connection  between  the  first  moment  of  the  popu¬ 
lation  distribution  ai  he  experimental  observable.  Thus, 
the  experimental  approach  can  be  based  on  measuring  prop¬ 
erties  that  depend  on  the  individual  moments  of  the  distribu¬ 
tion  function  sought. 

Some  experimental  measures  that  have  not  been  consid¬ 
ered  here  are  the  chemical  reaction  branching  ratios  to  form 


-.trengthened,  resulting  in  changes  in  vibrational  frequen-  the  isomeric  ethylene  products.  Branching  ratio  measure- 
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FIG  7.  Calculated  values  for  (£ '  as  a  function  of  (n)  Note  that  the  average 
thermal  vibrational  energy  of  TFE  at  300  K  is  about  470  cm 

ments  are  quite  feasible  in  many  systems  and,  in  many  cases, 
can  give  useful  information.  In  this  system,  however,  the 
reaction  threshold  energies  are  so  close  together  that  the  rel¬ 
ative  yields  are  insensitive  to  small  variations  of  the  param¬ 
eters  in  Table  II.  For  example,  variation  of  the  parameter 
'Wfromm  =  15  tom  =  25  is  calculated  to  change  the  rela¬ 
tive  yield  of  1,1 -difluoroethylene  from  :=0. 14  to  ~0. 16.  This 
variation  is  too  small  to  be  useful  in  refining  the  parameter 
values  in  the  present  investigation. 

In  the  present  study,  the  strategy  of  using  a  combination 
of  experimental  techniques  to  measure  different  attributes 
has  resulted  in  information  about  the  population  distribu¬ 
tions  formed  by  IRMPA.  This  information  will  be  used  in 
subsequent  studies  of  collisional  energy  transfer 1  and  btmo- 
lecular  reactions.  In  the  following  paper,’  the  calculated  po¬ 
pulation  distributions  are  used  to  predict  infrared  emission 
intensities,  which  arc  then  compared  to  those  observed  in 
expenmen  ts. 
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TABLE  I:  Frequencies  tmd  Aatomomkiriea  (i»  ci'1) 

A.  Water'1 
Frequencies  («,) 


365663 

1594  59 

3755  79 

Anharmonicity  Matrix  (A„) 

-42.576 

-15.933 

-165  824 

-15.933 

-16.813 

-20  332 

-165.824 

-20  332 

-47  566 

B 

Formaldehyde 

i: 

Frequencies  (w,) 


These  expressions  define  the  integration  region  boundary,  and 
therefore  the  sampling  domain  completely  encloses  the  integration 
region  when  3  >  1  When  3  *  1 .  the  sampling  domain  and  the 
integration  region  are  identical,  producing  100CTc  sampling  effi¬ 
ciency. 

For  degrees  of  freedom  other  than  vibrations,  a  sampling 
boundary  that  closely  resembles  the  integration  region  boundary 
may  be  chosen  for  ease  of  calculation.  It  is  not  always  desirable 
to  include  high-order  effects  in  the  sampling  boundaries,  because 
of  the  computational  labor  involved.  When  carrying  out  the 
sampling,  we  will  include  the  higher  order  effects  in  calculating 
Jig)  (see  eq  7).  In  this  way,  couplings  among  vibrations,  vibra¬ 
tions/rotations.  Coriolis  coupling,  Fermi  interactions,  etc.,  are 
incorporated  in  the  Monte  Carlo  integration. 

Results  and  Discussion 

Monte  Carlo  Integration.  The  Monte  Carlo  sampling  algorithm 
and  sampling  domain  described  above  were  tested  for  accuracy 
by  comparisons  with  exact  counts  of  states  based  on  spectroscopic 
vibrational  constants  for  water  and  formaldehyde  presented  in 
Table  I.  Although  the  spectroscopic  data  were  obtained  for 
energies  only  up  to  about  10000  cm"1  and  they  may  not  be  ap¬ 
propriate  for  higher  energies,  we  have  used  them  at  energies  up 
to  40000  cm'1  solely  for  the  purpose  of  testing  the  Monte  Carlo 
technique.  The  exact  counts  were  performed  by  using  the 
Stein-Rabinovitch  method  with  a  grain  size  of  2.7  cm"1 

The  harmonic  oscillator  calculations  were  performed  by  using 
the  tabulated  frequencies  but  setting  all  of  the  anharmomcities 
equal  to  zero:  the  tests  using  anharmonic  separable  oscillators  were 
carried  out  by  keeping  the  diagonal  anharmomcities  but  neglecting 
the  off-diagonal  anharmonicities.  We  chose  not  to  employ  more 
accurate  local  mode  descriptions  of  the  C-H  stretching  vibra¬ 
tions,1*  because  such  descriptions  require  additional  calculations 
and  knowledge  of  the  off-diagonal  couplings  between  the  stretching 
modes  and  they  neglect  the  other  off-diagonal  couplings 

In  Figure  2.  the  test  results  are  presented  as  ratios  of  the  Monte 
Carlo  results  divided  by  the  exact-count  results.  The  error  bars 
were  calculated  according  to  eq  8.  and  they  represent  ±  i  u  In 
all  four  of  the  test  cases,  the  ratio  of  Monte  Carlo/  exact  is  equal 
to  unity  within  I  a  for  most  of  the  data  points  Fewer  than 
one-third  of  the  test  points  (  E  >  0)  deviate  more  than  ±  1  cr  from 
unity,  as  expected  on  the  basis  of  the  normal  error  distribution. 
In  other  tests  mot  shown),  the  results  were  equally  satisfactory, 
and  it  was  found,  as  expected,  that  the  standard  deviation  is 
inversely  proportional  to  the  square  root  of  the  number  of  trials 

Tests  were  also  performed  to  determine  whether  the  sampling 
domain  completely  enclosed  the  integration  volume  and  whether 
the  achievable  precision  was  satisfactory  The  sampling  boundary 
was  defined  as  described  above,  and  the  parameter  3  was  varied 
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FjraJji  S.h  1981 


t!  i.M 

a 


3.  100 

b 


2811  42 

1755.858  1500.32 

1170.224  5 

1861.30 

1250.565 

-28.95 

Anharmonicity  Matrix  ( XtJ ) 
1.15  -23.03  -10099 

-193  32 

-49  78 

0.90 

1  15 

-9  926  -8  26 

-7  199 

-17  23 

6  581 

-23  03 

-8.26  -1.64 

-1.769 

600 

-29  861 

-10.099 

-7  199  -1  769 

-3  157 

-13  35 

-2  36 

-193  32 
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-17  63 
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Figure  2.  Nccuracy  tests  ul  waicr.  harmonic  oscillator  model,  tbi 
water  separable  anharmonic  oscillator  model,  to  formaldehyde  har¬ 
monic  oscillator  model,  id)  formaldehyde  veparable  anharmonic  oscil¬ 
lator  model 
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Multidimensional  Monte  Carlo  Integration 


TABLE  II:  D epessdessce  o>  d  »nd  Number  of  Trial* 


energy 

J 

sum 

trials 

efficiency,  % 

40000 

1  00 

Water  (Coupled) 
(1215  ±  0.006)  x  I0J 

50000 

100 

40000 

1  01 

(1.234  ±  0.012)  x  |0J 

17  532 

91 

40000 

1-25 

(1  226  ±  0.019)  x  I03 

25000 

19 

10000 

1.01 

Formaldehyde  (Coupled) 

(6  320  ±  0.075)  x  10!  100000 

92 

20000 

1  01 

(1  573  ±  0.025)  x  104 

100000 

36 

30  000 

1  01 

(  1  392  ±  0.025)  x  105 

100000 

78 

40000 

1  01 

(8  361  ±  0.157)  x  10s 

100000 

68 

to  ensure  that  the  integration  volume  was  completely  enclosed  by 
the  sampling  domain  (Table  II).  Inspection  of  Table  II  shows 
that  for  3  >  I  0  the  sampling  domain  appears  to  completely  enclose 
the  integration  volume,  even  when  the  vibrational  energy  is  40000 
cm'1,  nearly  10000  cm"1  greater  than  the  dissociation  energy. 
W  hen  3  =  10.  the  Monte  Carlo  sampling  is  100^  efficient:  every 
sample  point  lies  within  the  integration  region.  When  3  >  1.0. 
the  efficiency  is  lower,  but  convergence  of  the  integral  is  assured, 
since  points  that  fall  on  the  boundary  are  guaranteed  to  be  counted 
correctly.  For  most  of  the  calculations.  3  was  set  equal  to  I  01. 
which  assures  convergence  of  the  integral  without  much  sacrifice 
of  efficiency. 

Table  II  also  indicates  the  level  of  uncertainty  that  is  achievable 
with  this  simple  technique  for  a  given  number  of  trials.  For  the 
three-dimensional  (3-D)  integration  carried  out  for  the  water 
molecule  at  40000  cm'1,  a  \%  relative  standard  deviation  required 
17  532  trials.  If  200000  trials  had  been  used,  the  relative  standard 
deviation  would  have  been  ~0.3%.  This  result  can  be  compared 
roughly  to  that  obtained  by  Farantos  et  al..4  who  used  stratified 
sampling  in  the  3-D  classical  integration  for  ozone  with  9800  cm'1 
of  vibrational  energy;  they  achieved  ~0.3%  relative  standard 
deviation  by  calculating  '200000  trials  and  using  more  than  150 
"strata".  This  rough  comparison  shows  that  the  present  simple 
method  gives  results  of  high  precision,  without  the  necessity  for 
stratified  sampling.  If  the  present  simple  method  were  to  be 
employed  with  a  stratified  sampling  scheme,  even  greater  precision 
could  be  achieved  with  a  fixed  number  of  trials. 

Sonseparable  Anharmonic  Oscillators  These  calculations  were 
carried  out  by  using  all  of  the  anharmonicities.  The  sums  of 
vibrational  states  for  nonseparable  H-O11  and  CH.01J  are  sig¬ 
nificantly  larger  than  the  corresponding  harmonic  and  anharmonic 
separable  molecules,  as  shown  in  Figures  3  and  4.  Although  each 
of  these  molecules  shows  significant  “local  mode"  character,  the 
couplings  among  all  of  the  vibrational  modes  are  quite  significant. 
Couplings  between  the  C-H  stretching  modes  can  be  incorporated 
into  a  local  mode  description :3  that  is  more  accurate  than  the 
“separable  model"  employed  in  this  paper,  but  the  local  mode 
model  still  neglects  couplings  between  bends  and  stretches,  w  hich 
are  of  fundamental  importance  in.  for  example,  isomerization 
reactions. 

It  is  clear  from  the  comparisons  in  Figure  3  that  harmonic 
models  give  good  approximations  (factors  of  2)  to  the  full, 
nonseparable  descriptions,  even  at  high  energy  and  for  strongly 
anharmonic  and  coupled  molecules.  Thus,  for  fitting  RRKM 
transition  states  (with  unknown  frequencies!  to  experimental  data 
for  a  restricted  temperature  or  energy  range,  little  is  gamed  by 
employing  the  coupled  model,  unless  the  potential  energy  surface 
is  known  Otherwise,  the  unknown  harmonic  freouencies  can  be 
"adjusted"  to  give  agreement  with  the  data.  Furthermore,  the 
computational  speed  of  the  Monte  Carlo  technique  is  far  less  than 
that  of  the  Beyer-Swinhart  and  Stein-Rabinovitch  alrogithms 

In  this  paper,  we  have  considered  only  discrete  vibrational  states, 
based  on  spectroscopic  constants.  For  this  purpose,  eq  14  describes 
the  sampling  boundary  For  other  degrees  of  freedom  and  other 
intermode  coupling  schemes,  eq  14  must  be  replaced  by  other 
expressions,  which  may  make  evaluation  of  the  sampling  boundary 
more  difficult  In  principle,  the  methods  described  here  can  be 
extended  lo  any  arbitrary  molecular  model,  including  classical 
descriptions  Several  groups  have  employed  anharmonic  coupling 
schemes  m  classical  models  of  molecules  and  transition  states  4 
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Figure  4.  Sums  of  states  for  nonseparable  anharmonic  oscillator  models. 

In  future  work,  such  models  will  be  treated  with  the  Monte  Carlo 
method  described  here. 

For  theoretical  calculations  of  rate  constants  from  potential 
energy  surfaces,  where  the  off-diagonal  couplings  are  known,  the 
Monte  Carlo  techniques  gives  highly  accurate  sums  of  states, 
regardless  of  the  type  of  degree  of  freedom  or  how  a  is  coupled 
to  other  motions  of  the  molecule.  Harmonic  models  are  simply 
not  adequate,  since  they  give  sums  of  states  that  differ  significantly 
from  the  nonseparable  model.  The  mistaken  use  of  harmonic 
models  could  lead  to  significant  misinterpretations,  if  tunneling, 
avoided  crossings,  or  some  other  nonclassical  process  is  invoked 

In  general,  open  reaction  channels  introduce  dramatic  anhar¬ 
monic  coupling,  which  is  usually  nonseparable;  accurate  calcu¬ 
lations  require  nonseparable  models  under  these  circumstances 

1 1*1  For  example.  >ce  la)  Bhuiyan.  L.  B  .  Hase,  W  L.  J  Otrm  Ph\  t 
14*3  <0<2  Ibl  Wardlaw  D  M  Marcus.  R  X  /  Otem  Ptns  1*85. 

O.  Uis;  and  references  cued  therein 
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Thus,  the  Monte  Carlo  techniques  may  find  wide  use  in  describing 
van  der  Waals  complexes,  energy-transfer  'complexes',  and 
transition-state  properties  at  and  above  the  critical  energies  for 
reaction,  especially  when  the  reaction  is  studied  over  w  ide  energy 
and  temperature  ranges. 

Conclusions 

A  general  Monte  Carlo  technique  for  high-dimensional  integrals 
is  described  and  applied  to  calculations  of  sums  of  states  for 
nonseparable  anharmomc  molecules.  The  method  can  be  applied 
to  sampling  domains  with  complicated  boundaries,  and  it  has  an 
efficiency  that  can  equal  unity  in  realistic  applications  When 
applied  to  the  spectroscopically  determined  properties  of  real 
molecules,  the  results  show  that  the  calculational  technique  gives 
accurate  results  and  relative  standard  deviations  comparable  with 


tho,e  obtained  by  a  modern  stratified  sampling  method  which 
is  much  more  complicated  to  implement. 

The  sums  of  states  for  the  nonseparable  v,brations  of  H.O  and 
CH;0  are  significantly  different  from  those  for  the  separable 
vibrations  Harmonic  vibrations  are  probably  acceptable  for 
empirically  fitting  RRKM  models  to  experimental  data  but 
theoretical  calculations  of  rate  constants  and  lifetimes  must 
properly  account  for  the  nonseparable  characteristics  of  the  po¬ 
tential  energy  surface  The  Monte  Carlo  technique  described  here 
is  suitable  for  such  calculations. 
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Infrared  emission  spectra  were  obtained  for  1.  I,  2-tnfluoroethane  (TFE)  excited  by  infrared 
multiphoton  absorption  ( 1079.85  cm  ~  '|.  The  emission  features  show  that  the  HF  reaction 
product  is  formed  in  vibrational  states  up  to  about  v  =  3.  Furthermore,  emission  attributed  to 
F-CsC-H  was  observed  near  3320  cm"1,  indicating  that  the  difluoroethylene  primary  products 
of  TFE  decomposition  undergo  secondary  photolysis;  since  the  difluoroethylene  products  at 
room  temperature  do  not  absorb  laser  light,  they  must  be  formed  vibrationally  excited.  The 
emission  from  the  C-H  stretch  modes  of  TFE  was  readily  identified  near  2980  cm  ~ '  and  the 
emission  intensity  was  obtained  as  a  function  of  laser  fluence.  These  data  are  in  excellent 
agreement  with  predictions  based  on  the  theoretical  expression  for  fluorescence  intensity  and  the 
reconstructed  populations  determined  by  the  Master  Equation  calculations  described  in  the 
preceding  paper.  These  results  provide  additional  support  for  the  accuracy  of  the  reconstructed 
population  distributions  and  for  the  theory  relating  infrared  fluorescence  intensity  to  total 
vibrational  energy  in  polyatomic  molecules. 


INTRODUCTION 

In  the  preceding  paper,1  results  from  two  experimental 
techniques  were  used  in  conjunction  with  a  Master  Equation 
model2  to  determine  the  population  distributions  produced 
by  infrared  multiphoton  absorption  |  MPA |.  The  experimen¬ 
tal  techniques  used  were  optoacoustic  measurement  of  ab¬ 
sorbed  laser  energy2  and  VLP<t>4  measurement  of  infrared 
multiphoton  decomposition  tIRMPDl  yields.5  Although 
these  techniques  enabled  us  to  reconstruct  the  apporoximate 
population  distributions  produced  by  MPA.  they  do  not 
probe  directly  the  population  of  molecules  excited  to  ener¬ 
gies  below  the  decomposition  threshold. 

In  situ  experimental  techniques  that  directly  monitor 
the  vibrational  energies  of  highly  vibrationally  excited  mole¬ 
cules  in  their  electronic  ground  states  are  relatively  scarce. 
For  small  molecules,  laser  induced  fluorescence*  and  Ra¬ 
man  spectroscopic  methods'  are  appropriate,  but  for  large 
molecules  with  high  densities  of  states.  UV  absorbance""" 
and  infrared  fluorescence12"14  iIRFl  have  proved  to  be  effec¬ 
tive  Each  of  these  methods  requires  a  calibration  curve  that 
relates  the  observed  signal  to  the  vibrational  energy 

For  the  IRF  techniques,  the  calibration  curve  relating 
intensity  to  vibrational  energy  can  be  calculated  from  v  ibra- 
tional  spectroscopic  theory,  or  it  can  be  derived  from  expen- 
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ments.  An  important  motivation  for  the  present  experiments 
is  to  test  the  theoretical  predictions  of  IRF  intensity  against 
measurements  performed  using  MPA  to  produce  known 
populations  of  excited  molecules.  Validation  of  the  theoreti¬ 
cal  expression  in  the  present  experiments  supports  theoreti¬ 
cal  calculations  of  infrared  emission  spectra  of  highly  excited 
molecules.  Such  calculations  find  many  uses,  including  sim¬ 
ulations  of  IRF  from  interstellar  molecules  excited  by  ab¬ 
sorption  of  starlight.15 

The  experimental  approach  in  the  present  work  is  to 
obtain  time-  and  wavelength-resolved  IFF  data  as  functions 
of  laser  fluence  and  collider-gas  pressure.  By  extrapolating 
the  decay  curves  to  the  time  origin,  the  initial  IRF  intensity 
can  be  compared  to  that  predicted  by  vibrational  spectro¬ 
scopic  theory  and  the  Master  Equation  described  in  paper  I. 
In  the  third  paper  of  this  senes,1*  the  calibration  curve  is 
used  to  monitor  collisional  energy  transfer  and  determine 
the  energy  and  temperature  dependences  of  the  energy  trans¬ 
fer  step  size. 

EXPERIMENTAL 

Most  of  the  experimental  details  are  desenbed  in  paper 
I  Briefly,  high  power  TEA  CO-  lasers  operating  on  the  9  6 
nm  R  i22l  line  1 1079  85  cm  "  'l  were  used  to  excite  1.  1.2- 
tnfluoroethane.  Laser  energy  was  measured  with  a  volume- 
absorbing  calonmeter  and  fluence  was  calculated  from  the 
measured  energy  and  the  laser  beam  cross-sectional  area 
Gas  pressures  were  measured  with  a  capacitance  manome¬ 
ter 
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The  cell  used  for  the  IRF  experiments  was  the  same  as 
that  used  for  the  optoacoustic  measurements  in  paper  I.  The 
fluorescence  was  viewed  through  a  KC1  side  window  by  a 
liquid  nitrogen  cooled  two-color  infrared  detector  (Infrared 
Associates,  Inc.  I,  which  consists  of  an  InSb  detector  overly¬ 
ing  a  HgCdTe  photodiode.  Each  detector  layer  is  equipped 
with  a  separate  preamplifier  and  biasing  circuit.  The  InSb 
layer  is  sensitive  from  1  to  5.5/im;  at  longer  wavelengths  it  is 
transparent,  allowing  the  HgCdTe  detector  to  respond  to 
wavelengths  from  5.5  to  16/xm.  In  the  present  experiments, 
only  the  InSb  portion  of  the  detector  was  used.  The  original 
time  response  of  the  InSb  detector/preamplifier  was  1-2  ps, 
but  it  was  susceptible  to  electrical  interference  from  the  la¬ 


sers;  thus  it  was  modified  to  give  a  slower  time  response  (5-fi 
/as)  and  less  susceptibility,  although  there  is  still  some  evi¬ 
dence  of  electrical  “ringing"  With  this  modification  and  due 
to  the  ringing,  accurate  measurements  were  limited  to  decay 
rate  constants  slower  than  about  7x  104  s~'.  The  signal 
from  the  InSb  detector/preamplifier  was  captured  using  a 
transient  recorder  iBiomation  805)  and  signal  averager  iNi- 
colet  i072)  interfaced  to  a  laboratory  computer  isee  paper  I 
for  details). 

Wavelength  resolution  was  achieved  using  a  three-seg¬ 
ment  circular  variable  filter  lOCLII  mounted  in  a  custom 
Dewar  with  the  detector  and  cooled  to  77  K.  The  slits  were 
1 . 5  mm,  giving  a  resolution  of  about  2%  of  the  center  wave- 


FIG  2  Infrared  fluorescence  decay  curve 
20  mTorr  TFE  and  20  mTorr  argon 
'P  -  I  HO  J  cm  253  laser  shots.  3frt>5 
cm 
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FIG  3  Infrared  fluorescence  decay  curve. 
19  mTorr  TFE  and  19  mTorr  argon. 
4>  =  1.95  J  cm  121  laser  show,  3020 

cm"1. 


length.  This  resolution  is  sufficient  to  distinguish  C-H 
stretch  emission  bands  from  HF  emission  and  other  interfer¬ 
ences.  The  filter  bandpass  and  relative  wavelength  calibra¬ 
tion  i  wavelength  vs  rotation  anglel  were  determined  using  an 
infrared  spectrometer.  Absolute  positioning  of  the  rotation 
angle  was  achieved  by  observing  the  CO;  laser  line  at  9  26 
urn  and  fluorescence  from  COz  near  4.2  pm.  We  estimate  the 
absolute  wavelength  calibration  to  be  accurate  to  within 
-  20  cm  "  '. 

RESULTS  ANO  DISCUSSION 

A  surv  ey  of  the  IRF  emission  was  earned  out  for  wave¬ 
lengths  from  about  2000  to  4000  cm-'  and  laser  fluences 


from0.6to  1.95  Jem-2.  In  these  experiments,  a  1/1  mixture 
of  1 , 1 , 2-tnfluoroethane  (TFE)  in  argon  at  a  total  pressure  of 
32-40  mTorr  flowed  slowly  through  the  cell.  Argon  was 
added  to  slow  the  rates  of  diffusion  and  thermal  conductiv¬ 
ity,  and  to  moderate  the  temperature  excursions.  The  laser 
pulse  repetition  frequency  was  0.2-8  Hz,  depending  on  laser 
fluence,  and  care  was  taken  not  deplete  the  TFE  by  more 
than  a  few  percent.  Depending  on  signal  strength,  between 
50  and  1 5  000  laser  shots  were  accumulateed  for  each  IRF 
decay  curve  to  improve  the  signal/ noise  ratio. 

Fluorescence  decay  curves  were  obtained  at  wavelength 
intervals  of  50-100  cm  - 1  and  a  few  representative  examples 
are  presented  in  Fig.  1-3.  Figure  3  shows  a  relatively  simple 
exponential  decay,  while  Figs.  1  and  2  show  more  complicat- 


FIG  4  Time-resolved  infrared  emission 
■.pectrum  20  mTorr  TFE  and  20  mTorr 
argon,  =  0  59  J  cm  ' :.  estimated  uncer¬ 
tainties  are  -  I0T>  The  top  solid  line  is 
for  3  us  following  the  laser  pulse  and  each 
successive  line  is  for  additional  I0,us  delay 
times 
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FIG.  5.  Time- resolved  infrared  emission 
spectrum.  20  mTorr  TFE  and  20  mTorr 
argon,  ♦  =  1 . 1 7  J  cm  ’ 1,  estimated  uncer¬ 
tainties  are  ±  10%.  Note  the  different 
wavelength  scale  and  see  Fig.  4  for  key. 


ed  behavior  that  is  evidence  for  a  collisional  cascade.  All  of 
the  decay  curves  were  individually  least-squares  fitted  by 
sums  of  exponentials.  By  combining  the  fitted  results,  we 
obtained  the  time-  and  wavelength-resolved  spectra  present¬ 
ed  in  Figs.  4-6,  which  correspond  to  low,  medium,  and  high 
laser  fluences,  respectively. 

Inspection  of  Figs.  4-6  reveals  three  distinct  spectro¬ 
scopic  features  at  2980,  3320,  and  3500-3800  cm” in  addi¬ 
tion  Fig.  5  shows  continuum  emission  between  2100  and 


creases  dramatically  with  laser  fluence.  The  wavelength  cor¬ 
responds  to  that  of  the  P  branch  of  HF,  which  is  a  primary 
reaction  product18;  the  wavelength  range  of  the  CVF  is  not 
sufficient  to  observe  the  Q  branch.  This  feature  changes 
width  and  center  frequency  as  a  function  of  time,  indicating 
that  several  vibrational  levels  of  HF  are  initially  formed  in 
the  decomposition  of  TFE  and  a  vibrational  cascade  ensues. 
The  detailed  shapes  of  the  decay  curves  in  Figs.  1  and  2  are 
explained  by  the  shift  and  broadening  of  the  3800  cm  ” 1  fea- 


Zellweger.  Brown,  and  Barker-  Vibrationally  excited  populations.  II 


6265 


the  decomposition  reaction.  A  more  detailed  analysis  is  be¬ 
yond  'he  scope  of  the  present  work. 

The  origin  of  the  feature  at  3320  cm" 1  is  uncertain,  but 
it  only  appears  at  very  high  fluences  and  is  about  the  same 
frequency  as  the  C-H  stretching  mode  in  acetylene. ;o  Thus, 
we  tentatively  conclude  that  it  arises  from  vibrationally  ex¬ 
cited  F-C==C-H  produced  by  secondary  photolysis  of  di- 
fluorethylene  products  from  the  primary  dissociation: 

CF.H-CFH,  -  nhv—  CF;  =  CH;  -r  HF  iA) 

—cis-  or  rra/is-CFH  =  CFH  •+■  HF  t  B) 

—  CF-CFH,  +  HF 

—cis-  or  ,'raru-CFH  =  CFH  -t-  HF  ,  |C) 

Difluoroethylene|vib.|  mhv— -F-C=C-H  -f  HF  .  iD) 

The  secondary  photolysis  [process  |D)]  can  only  occur 
with  vibrationally  excited  molecules,  because  the  room  tem¬ 
perature  difluoroethy lene  products  do  not  absorb  at  the  laser 
wavelength,  as  discussed  in  paper  I.  It  is  interesting  to  note, 
then,  that  both  products  of  the  photolysis  must  contain  sig¬ 
nificant  amounts  of  vibrational  energy.  Because  of  the 
strained  three-  and  four-center  cyclic  transition  state  struc¬ 
tures,  it  is  expected  that  vibrationally  excited  HF  and  di- 
fluoroethylene  will  be  produced,  but  the  extent  of  excitation 
and  the  apportionment  of  energy  varies,  depending  on  the 
particular  reaction  channel.18  For  example,  the  slightly 
higher  energy  three-center  pathway  (C)  requires  isomeriza¬ 
tion  of  the  carbene  primary  product  to  form  difluoroethy- 
lene  and  the  difluoroethylene  produced  by  this  channel  is 
apportioned  more  vibrational  excitation  than  the  four-cen¬ 
ter  channels. 

The  spectroscopic  feature  at  2980  cm  “ 1  corresponds  to 
the  C-H  stretching  modes  in  TFE.21  The  C-H  stretching 
modes  in  cis-  and  rra/w-difluoroethylene  have  frequencies20 
near  3100  cm-1  and  cannot  contribute  to  the  2980  cm-1 
feature.  The  monotonic  decay  curves  at  this  frequency  (see 
Fig.  3 1  indicate  relatively  simple  collisional  deactivation  of 


the  excited  TFE  and  they  can  be  accurately  fitted  with  a 
single  exponential  decay.  The  empirical  fitting  on  the  decay 
curves  is  used  in  the  present  work  only  for  the  purpose  of 
extrapolation  to  obtam  the  initial  intensity;  collisional  deac¬ 
tivation  of  TFE  is  discussed  in  paper  III. 16  The  experimental 
IRF  decay  curves  for  2980  cm  ~ 1  were  fitted  by  a  single  expo¬ 
nential  with  the  Bevington  version  of  the  Marquardt  nonlin¬ 
ear  least  squares  algorithm.22  The  first  6  ps  of  the  decay 
following  the  laser  pulse  was  neglected  (see  above)  and  the 
least-squares  result  was  extrapolated  to  r  =  0  in  order  to  ob¬ 
tain  the  initial  IRF  intensities.  The  initial  intensities  were 
normalized  according  to  the  partial  pressure  of  TFE,  and  the 
results  are  presented  vs  laser  fluence  in  Fig.  7. 

For  comparison  with  the  experimental  results,  the  IRF 
intensity  was  calculated  by  using  the  population  distribu¬ 
tions  reconstructed  in  paper  I  and  the  theoretical  expression 
for  the  microcanonical  IRF  intensity.  The  theoretical 
expression  relating  IRF  intensity  to  vibrational  energy  m 
polyatomic  molecules  has  been  discussed  elsewhere  in  de¬ 
tail.21213  The  expression  is  written 

Mode* 

,{E )  =  .V,  X  £  hv,  Aj'0  £  vp,  vhv,)/p,\E),  (1) 

j  -  I  1 

where  S,  is  the  number  of  excited  molecules,  A  1,0  is  the 
Einstein  coefficient  for  emission  from  the  fundamental  of  the 
ith  mode,  u  is  the  occupation  number  of  the  fth  mode,  vmUL  is 
the  maximum  occupation  number  for  vibrational  energy 
E,  p,  [E )  is  the  density  of  vibrational  states  for  all  s  vibration¬ 
al  modes,  and  p,  _  ,  \E  —  vhv, )  is  the  density  of  states  for  the 
s  -  1  modes,  not  including  the  emitting  mode  and  the  energy 
contained  in  it. 

Briefly,  the  theoretical  expression  is  based  on  the  funda¬ 
mental  relationship  between  IRF  intensity  and  vibrational 
quantum  number,23  and  it  incorporates  two  assumptions:  ( 1 1 
the  harmonic  oscillator  approximation  relating  Einstein  co¬ 
efficients  for  higher  vibrational  levels  to  that  for  the  funda¬ 
mental  and  1 2)  the  ergodic  assumption,  which  states  that  the 
vibrational  energy  is  statistically  distributed  among  the  vi- 
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from  C-H  modes  to  vibrational  energy  Thu 
curve  was  calculated  using  theoretical  Eq.  1 1 1, 
see  the  text  for  details. 


brational  modes  on  the  time  scale  of  the  fluorescence  i  micro¬ 
seconds).  The  accuracy  of  this  expression  also  relies  on  the 
densities  of  states  calculations,  which  can  be  earned  out  us¬ 
ing  exact-count  algorithms. 

For  the  theoretical  prediction  of  IRF  intensity  vs  vibra¬ 
tional  energy,  the  Stein-Rabinovitch  method  of  exact 
counts24  was  used  to  obtain  the  densities  of  states,  based  on 
the  vibrational  assignments  for  TFE.21  The  hindered  inter¬ 
nal  rotation  was  treated  as  a  vibration  in  these  calculations. 
The  Einstein  coefficients  for  the  fundamental  transitions  are 
not  known  and  they  are  included  in  the  undetermined  scale 
factor,  which  also  depends  on  the  experiment  geometrical 
configuration,  detector  sensitivity,  amplifier  settings,  etc.  In 
Fig.  8,  the  predicted  microcanonical  dependence  of  IRF  in¬ 
tensity  on  vibrational  energy  is  presented  in  arbitrary  units. 

To  compare  the  predicted  IRF  intensities  with  those 
observed  at  different  laser  fluences,  the  population  distribu¬ 
tions  from  paper  I  were  used  with  the  microcanonical  IRF 
intensities  to  obtain  ensemble-averaged  IRF  intensities  vs 
laser  fluence.  These  theoretical  results  are  shown  in  Fig.  7  as 
the  dashed  line,  which  has  been  scaled  to  match  the  magni¬ 
tude  of  the  experimental  data  at  a  fluence  of  1 .0  J  cm  ~ 2.  The 
agreement  between  experiment  and  theory  is  excellent. 

The  excellence  of  the  agreement  between  experiment 
and  theory  may  be  viewed  in  several  ways.  If  Eq.  (II  is  ac¬ 
cepted  as  accurate,  the  agreement  provides  additional  sup¬ 
port  for  the  reconstructed  population  distributions  obtained 
in  paper  I.  This  is  a  legitimate  point  of  view  in  that  Eq.  1 1 )  has 
been  successfully  tested  is  several  other  systems  on  other 
molecules,1211  supporting  the  inference  that  it  should  accur¬ 
ately  describe  IRF  in  the  molecule  TFE.  The  other  point  of 
view  is  that  Eq.  ( 1)  must  be  tested  by  experiment,  despite  the 
elegance  of  its  derivation.  The  reconstructed  population  dis¬ 
tributions  are  derived  from  experiments  in  paper  I  and  ther- 
fore  the  experiments  reported  here  provide  an  experimental 
test  and  demonstration  of  the  accuracy  of  Eq.  ill,  when  ap¬ 
plied  to  TFE. 

Either  of  these  points  of  view  leads  to  the  conclusion 


that  the  Master  Equation  from  paper  I  and  Eq.  ( 1 )  constitute 
a  self-consistent  description  of  the  results  from  three  inde¬ 
pendent  experimental  techniques.  Considering  the  highly 
constrained  parametrization  of  the  Master  Equation  as  de¬ 
scribed  in  paper  I  and  the  fact  that  Eq.  (1)  accurately  de¬ 
scribes  IRF  in  other  systems  as  well  as  in  this  one,  we  con¬ 
clude  that  the  reconstructed  population  distributions  are 
good  representations  of  the  actual  distributions  and  that 
I RF  can  be  confidently  used  to  monitor  vibrational  energy  in 
TFE  excited  by  multiphoton  absorption.  This  technique  is 
used  in  paper  III  to  determine  energy-transfer  properties  of 
TFE. 
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Infrared  multiphoton  absorption  at  several  laser  fluences  (1079  85  cm'1)  has  been  used  to  generate  vibrationally  excited 
1 .  1 ,2-trifIuoroethane  molecule*  at  energies  below  the  reaction  threshold.  A  technique  is  demonstrated  in  which  collisional 
deactivation  of  the  vibrationally  excited  molecules  by  argon  was  monitored  with  time-resolved  infrared  fluorescence  from 
the  C-H  stretch  mode*.  The  experimental  results  were  simulated  by  using  a  full  collisional  master  equation  treatment  that 
incorporated  Quack's  theory  of  multiphoton  absorption,  as  parameterized  elsewhere  for  TFE.  and  a  calibration  cuive  relating 
fluorescence  to  vibrational  energy  (previously  determined  for  TFE).  When  the  exponential  model  for  collision  transition 
probabilities  was  used,  the  simulations  gave  ( ±E)t  =*  (200  ±  20)  +  (0  005  ±  0.002)£.  where  ( A£)d  is  the  energy  step  size 
(cm'1)  for  deactivating  collisions  by  argon  and  E  is  the  vibrational  energy  of  the  excited  TFE.  The  simulations  showed  that 
( SE)t  has  a  significant  energy  dependence,  but  it  does  not  exhibit  a  detectable  temperature  dependence  over  the  range  from 
about  400  to  1000  K.  In  the  introduction,  current  research  activities  by  several  workers  active  in  the  area  of  large  molecule 
energy  transfer  are  briefly  surveyed. 


IntroductiM 

Recently,  three  major  questions  and  several  lesser  ones  have 
arisen  about  vibrational  energy  transfer  involving  large  molecules. 
These  question*  are  all  concerned  with  ( A£>d.  the  average  amount 
of  energy  transferred  in  deactivating  collisions,  which  is  related 
by  microscopic  reversibility  and  detailed  balance  to  <  A£>,  the 
average  amount  of  energy  transferred  in  both  activating  and 
deactivating  collisions.1  The  major  questions  are  concerned  with 
( I )  the  magnitude  of  ( A£)d,  (2)  whether  (A£>„  depends  on  the 


’  All  of  the  experimental  work  wet  carried  out  SRI  International. 


vibrational  energy  of  the  highly  vibrationally  excited  molecule, 
and  (3)  how  ( A£)d  varies  with  the  collider  bath  gas  temperature. 
Few  energy-transfer  systems  have  been  examined  in  enough  detail 
to  address  any  of  these  questions  and  the  conflicting  results  may 
stem  from  the  different  molecules  studied,  or  from  undetected 
experimental  artifacts. 

One  of  the  best  approaches  to  determining  the  actual  reliability 
of  experimental  data  is  to  compare  results  from  several  different 


1 1 )  For  an  interesting  discussion  of  the  relative  menu  of  citing  ( A£>  or 
<±£)t.  see  Gilbert.  R.  G.  Chtm  Phyi  Lett  I9M.  S3.  21. 
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experimental  techniques.  This  was  one  of  the  primary  motivations 
of  our  work  on  azulene.2  which  employed  visible  and  UV  laser 
photoexcitation  to  prepare  the  vibrationally  excited  molecules  at 
selected  energies  and  time-resolved  infrared  fluorescence  to 
monitor  their  collisional  deactivation.  These  experiments  on  a 
nonreactive  system  were  intended  to  test  and  extend  the  large  body 
of  energy-transfer  data  obtained  in  unimolecular  reaction  rate 
studies.’ 

Recently,  Hippier.  Troe,  and  co-workers'*  have  also  investigated 
photoexcited  azulene  at  one  excitation  energy,  but  they  monitored 
the  deactivation  by  time-resolved  UV-absorbance  techniques  and 
reached  slightly  different  conclusions.  Specifically,  both  studies 
agree  that  the  magnitudes  of  <A£)d  in  the  azulene  system  are 
somewhat  smaller  for  efficient  colliders  than  estimated  on  the  basis 
of  unimolecular  reaction  rate  studies,  but  Hippier.  Troe,  and 
co-workers  concluded  that  (A£>d  is  essentially  independent  of 
vibrational  energy  content,  contrary  to  our  conclusions,  but  in 
agreement  with  their  own  work  on  other  molecules.’ 

Not  all  UV-absorbance  experiments  give  energy-independent 
values  for  (A £)d.  Experiments  by  Troe  and  co-workers6  on 
tnatorruc  molecules  give  values  for  ( A£)d  that  depend  dramatically 
on  internal  energy.  Also,  in  time-resolved  UV-absorbance  ex¬ 
periments  carried  out  by  Nakashima  and  Yoshihara  on  benzene7 
and  hexafluarobenzene,*  the  results  indicate  that  <  A£)d  is  directly 
proportional  to  internal  energy.  Thus  experiments  using  infrared 
fluorescence  on  these  molecules  would  provide  another  useful 
comparison. 

In  a  second  controversy  concerning  the  properties  of  (A£)d, 
different  studies  have  arrived  at  significantly  different  conclusions 
about  its  temperature  dependence.  It  was  concluded  on  the  basis 
of  infrared  fluorescence  experiments2*  that  < A£>d  shows  a  slight 
decrease  in  magnitude  as  the  collider  bath  gas  temperature  is 
raised.  On  the  other  hand,  multiple-channel  unimolecular  reaction 
studies  indicate  a  very  strong  decrease  in  magnitude  of  <A£)d 
with  temperature.*  while  UV-absorbance  experiments  indicate  a 
weak  temperature  dependence  that  can  change  sign  depending 
on  the  efficiency  of  the  collider.10  There  is  no  clear  explanation 
for  these  differences,  and  new  experimental  techniques  may  help 
to  resolve  the  conflicts. 

In  a  number  of  different  laboratories,  new  experimental  tech¬ 
niques  for  energy-transfer  studies  are  under  development.  For 
example,  photoexcitation  based  on  laser  pumping  of  high  overtones 
of  the  C-H  stretching  vibrations  in  potentially  reactive  molecules 
has  recently  been  applied  in  efforts  to  determine  energy-transfer 
properties."12  Also,  several  versions  of  crossed  molecular  beam 
experiments  are  under  study  in  at  least  three  different  labora¬ 
tories. 1  ^ 1  ’  Time-resolved  optoacousbc  methods  also  show  promise 

(2)  la)  Smith.  O  P ;  Barker.  J.  R.  Chem.  Phys.  Lett.  19*1,  78,  253.  (b) 
Roaai.  M  J  :  Barker.  J.  R.  Chem  Phys.  Lett.  19*2. 85.  21.  (c)  Barker.  J.  R ; 
Roui.  M  J  ;  Pladziewicz,  J.  R.  Chem.  Phys.  Lett.  IN},  90,  99.  (d)  Roaai, 

M  J  .  Pladziewicz.  J.  R.;  Barker,  J.  R.  J  Chem.  Phys.  19»3,  78,  6695  (e) 
Barker.  J  R  J  Phys.  Chem.  I9S4,  88.  II.  (f)  Barker.  J.  R  ;  Golden.  R.  E. 

J  Phys  Chem.  19*4. 88.  1012.  (|>  Font,  W.;  Barker.  1.  R.  J  Chem  Phvs 
19*5.  83,  124. 

(3)  (a)  Tardy.  D.  C..  Rabinovitch.  B.  S.  Chem.  Rev  19T7,  77,  369  (b) 
Quack.  M.;  Troe,  J.  Spec.  Period.  Rep.,  Chem.  Soc.  1977,  2,  175. 

(4)  (a)  Hippier,  H.  Her.  Bunsenges.  Phys.  Chem  19*5,  89.  304  (b) 

Hippier.  H.;  Lindemann.  L.;  Troe,  J  J.  Chem.  Phys  19*5.  83.  3906 

(5)  (a)  Hippier,  H.;  Troe,  J  .  Wendelken.  H  J  Chem  Phvs  Lett  19*1. 
84.  257  (b)  Hippier,  H  .  Troe.  I  .  Wendelken.  H.  J  J  Chem  Phvs  19*3, 

6709  (c)  Hippier,  H  ,  Troe,  3  .  Wendelken,  H.  J  J.  Chem.  Phys  19*3. 
'8.  6718. 

(6)  (a)  Dove,  J  E.;  Hippier,  H  .;  Troe.  J.  J  Chem  Phys.  19*5,  82.  1907 
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for  making  contributions  to  our  knowledge  of  energy  transfer 
involving  large  molecule*. 16  In  addition,  the  disposal  of  vibrational 
energy2*  in  the  energy-transfer  act  is  being  investigated  by  cold- 
gas-filter  experiments17  and  by  time-resolved  laser  diode  absorption 
measurements. 11  The  development  and  application  of  these 
techniques  promise  to  improve  vastly  knowledge  of  energy  transfer. 
In  the  following  sections,  we  will  describe  the  use  of  infrared 
multiphoton  absorption  (MPA)  to  prepare  ensembles  of  vibra¬ 
tionally  excited  molecules,  which  can  be  monitored  by  infrared 
fluorescence  (IRF). 

Multi  photo*  Absorption/Infrared  Fluorescence  Technique 

Infrared  multiphoton  absorption  is  fairly  general,  in  that  many 
molecules  absorb  light  from  high-power  C02  lasers.  A  large 
fraction  of  the  irradiated  molecules  are  excited  by  the  laser  light, 
but  the  resulting  distribution  function  is  broad  and  must  be 
characterized,  as  in  the  First  paper  of  this  series.'*  The  earliest 
applications  of  the  infrared  pumping  technique  analyzed  the  effects 
of  pressure  on  multiphoton  dissociation  yields.20-21  The  analysis 
is  quite  complex,  because  the  laser  pulse  shape  and  the  extreme 
temperature  excursions  after  the  pulse  both  play  important 
roles.21-22  Moreover,  the  interpretation  depends  critically  on 
accurate  knowledge  of  the  specific  unimolecular  rate  constants, 
which  must  be  determined  in  separate  experiments.  Nonetheless, 
the  energy-transfer  results  showed  considerable  promise. 

The  present  energy-transfer  experiments  differ  from  the  earlier 
dissociation  yield  techniques  in  several  important  ways.  The  most 
important  difference  is  that  the  energy-transfer  measurements 
do  not  depend  on  reaction  yield  measurements,  but  on  infrared 
fluorescence  from  nonreacting  molecules.  Also,  milder  conditions 
are  used,  leading  to  much  smaller  temperature  excursions,  easing 
interpretations  of  the  results.  There  is  the  added  benefit  that  the 
laser  pulse  shape  is  not  important  in  the  data  analysis  and  uni¬ 
molecular  rate  parameters  need  not  be  known  with  great  accuracy. 
However,  the  dependence  of  infrared  fluorescence  intensity  on 
internal  vibrational  energy  must  be  known  from  experiment  or 
from  theory.23  In  the  second  paper  of  this  series,  the  theoretical 
expression  relating  IRF  intensity  to  vibrational  energy  was  com¬ 
pared  to  experimental  results  obtained  from  infrared  multipboton 
absorption,  and  the  comparison  showed  very  good  agreement.24 
Thus,  the  basis  for  energy-transfer  experiments  has  been  estab¬ 
lished. 

The  approach  in  these  experiments  is  to  use  multiphoton  ab¬ 
sorption  to  prepare  nonreactive  ensembles  of  vibrationally  excited 
1,1.2-trifluoroethane  (TFE)  molecules  (the  initial  population 
distributions  are  known  with  good  accuracy,  based  on  the  results 
of  paper  1).  The  collisional  deactivation  of  the  excited  molecules 
is  monitored  by  time-resolved  IRF  from  the  C-H  stretching  modes. 
By  varying  the  gas  pressure,  the  deactivation  rate  can  be  varied 
systematically,  and  by  changing  gas  composition,  the  temperature 
excursions  can  also  be  controlled.  Thus,  investigation  of  the  effects 
of  laser  fluence  and  of  gas  composition  enables  us  to  determine 
the  dependence  of  ( A£)d  on  vibrational  energy  and  on  temper¬ 
ature. 

Because  of  the  temperature  excursions  and  because  the  initial 
population  distributions  are  relatively  broad,  a  master  equation 
treatment  is  used  in  data  analysis.  The  master  equation  imple¬ 
mentation  used  to  simulate  the  experimental  data  has  been  de- 
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Figure  1.  Reconstructed  population  distributions.  These  population 
distributions  are  from  the  master  equation  treatment,  as  parametertzed 
in  paper  I.  For  comparison,  the  300  K  thermal  distribution  (X  1 , 2)  is 
shown  in  the  lop  panel. 


scribed  in  detail  elsewhere.2122  Recent  improvements  described 
in  paper  I  consist  of  using  exact-count  densities  of  states,  rather 
than  the  Whitten-Rabinovitch  approximation,  and  incorporating 
Quack's  theory  of  infrared  multiphoton  absorption.25  Paramet¬ 
erization  of  the  theory  to  fit  infrared  multiphoton  absorption  by 
TFE  allowed  us  to  reconstruct  the  population  distribution  produced 
by  a  selected  laser  fluence.'9 

Although  the  experimental  population  distributions  are  not 
uniquely  determined  by  the  master  equation  simulation  approach, 
the  theoretical  reconstruction  is  a  good  approximation  and  is  used 
as  input  in  the  energy-transfer  simulations.  It  was  shown  in  paper 
(I  that  the  reconstructed  populations  could  be  coupled  with  the 
theoretical  expression  for  IRF  dependence  on  vibrational  energy 
to  accurately  predict  the  observed  IRF  intensity  produced  by  laser 
fluencea  up  to  about  2  J  cm'1.1* 

Examples  of  the  reconstructed  population  distributions  for  TFE 
are  presented  in  Figure  I.  These  populations  show  the  broad 
distributions  often  predicted  by  most  models  for  infrared  multi¬ 
photon  absorption,  but  there  are  additional  features  that  are  of 
interest.  At  low  fluences,  significant  numbers  of  molecules  remain 
in  the  initial  thermal  distribution  and  are  not  pumped  up  the 
ladder.  This  feature  is  a  direct  consequence  of  the  behavior  of 
Quack  case  C,25  which  dominates  in  the  sparse  density  of  states 
region  at  low  energy.  At  very  high  fluences.  the  population 
distributions  tend  to  ‘stack  up"  at  high  energies  and  become  much 
narrower.  This  effect  is  due  to  a  pumping  ‘bottleneck’  caused 
by  the  anharmonic  shift  of  the  pumped  vibrational  mode  (a  C-F 
stretch).1*  The  anharmonicity  is  large  for  energies  near  the  re¬ 
action  threshold,  because  the  electronic  hybridization  of  the  carbon 


(24)  Paper  II:  Zellweger.  J.-M.;  Brown.  T.  C.;  Barker.  J.  R.  J.  Chtm. 
Phys.  IMS.  8S,  <261:  this  paper  describes  the  infrared  fluorescence  observed. 

(23)  Quack.  M.  Bet.  Buiutnfts.  Phys  Chtm.  IMI.  85,  318  (for  a  treat¬ 
ment  that  considers  fluence-dependent  fractionation  of  pumped  and  non- 
pumped  moienilea,  see  fang,  J.  C  ,  Setter.  D.  W.;  Dancn.  W.  C.  J.  Am.  Chtm. 
Soc.  IMS.  104.  3440). 


:o  I - I - ! - ! - 1 - - 

2:00  2S00  2900  3300  3~03 


WAVELENGTH  'cm  ’1 

Figure  2.  Infrared  emission  spectra.  From  the  experiments  in  paper  II. 
these  time-  and  wavelength-resolved  curves  correspond  to  about  20  mtorr 
of  TFE  diluted  in  20  mtorr  of  Ar  and  pumped  at  1079.83  cm'1  with  a 
laser  fluence  of  l  .2  J  cm'2.  The  first  spectrum  (top  solid  line)  is  8  ns  after 
the  pulse,  and  each  succeeding  spectrum  occurs  10  ys  later. 


Figure  3.  Infrared  fluorescence  decay  curve.  Emission  at  2980  cm'1  from 
9  mtorr  of  TFE  diluted  in  81  mtorr  of  Ar  and  pumped  with  a  laser 
fluence  of  I  07  J  cm'2;  the  least-squares  decay  rate  constant  is  (4  09  ± 
0.22)  x  10*  s'1. 


atoms  changes  from  sp3  to  sp2  as  the  molecule  reacts  by  HF 
molecular  elimination. 

Experimental  Section 

The  experimental  apparatus  and  techniques  are  the  same  as 
those  used  in  papers  1  and  11.  Briefly,  a  high-power  pulsed  CO, 
laser  tuned  to  the  1079.85-cm'1  transition  is  used  to  excite 
1 . 1 ,2-triiluoroethane  (TFE)  diluted  in  argon.  Infrared  fluorescence 
from  the  C-H  stretch  modes  of  TFE  is  used  to  monitor  the 
cotlisional  deactivation  of  the  excited  molecules.  The  IRF  ts 
spectrally  isolated  with  a  circular  variable  Filter  set  at  2980  cm  1 
The  resolution  of  the  filter  at  this  wavelength  is  about  60  cm'1, 
sufficient  to  cleanly  distinguish  the  C-H  IRF  from  HF  fluores¬ 
cence  and  other  interferences,  as  shown  in  paper  11,  where  time- 
and  wavelength-resolved  IRF  spectra  were  obtained;  an  example 
is  presented  in  Figure  2.  The  IRF  is  monitored  by  using  the  InSb 
portion  of  a  two-color  infrared  detector  ( Infrared  Associates,  Inc.) 
and  the  signal  is  captured  with  a  transient  recorder /signal  averager 
interfaced  to  a  laboratory  computer.  The  overall  time  constant 
of  the  detector,  amplifiers,  and  digital  electronics  is  1-2  us.  but 
the  first  4-6  us  of  the  signal  is  obliterated  by  electrical  interference 
from  the  pulsed  laser. 

Results  sod  DUcasaioa 

A  typical  IRF  decay  is  presented  in  Figure  3,  which  shows  the 
effect  of  the  electrical  interference.  This  decay  curve  was  obtained 
with  1030  laser  shots  and  it  illustrates  that  the  signals  are  quite 
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LASER  -ljENCE  Ucm-*; 

Figure  4.  Infrared  fluorescence  decay  rate  constants.  The  four  solid 
curves  are  calculated  by  using  the  master  equation  for  9  mtorr  of  TFE 
diluted  in  argoa:  (O)  358.  (V)  225,  (A)  120.  and  (0)  81  mtorr  of  argon, 
from  top  to  bottom.  The  experimental  points  show  one  standard  devia¬ 
tion  uncertainties. 
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FINAL  TEMPERATURE  v«.  FLUENCE 

Figwt  S.  Final  temperature  after  equilibration.  Calculated  by  using  the 
master  equation  for  9  mtorr  of  TFE  diluted  in  argon.  From  top  to 
bottom,  the  four  curves  correspond  to  81.  120,  225,  and  358  mtorr  of 
argon. 

noisy,  even  with  extensive  signal  averaging.  Nevertheless,  a 
nonlinear  least-squares24  fit  of  the  data  beginning  6.6  ms  after  the 
laser  pulse  (which  occurs  at  29  ms)  gives  the  decay  constant  with 
an  uncertainty  (±l<r)  of  only  about  5%.  This  level  of  precision 
is  adequate  for  present  purposes,  but  it  is  not  sufficient  to  detect 
more  subtle  effects,  such  as  deviations  from  exact  exponential 
decay. 

To  demonstrate  the  application  of  the  technique,  a  series  of 
experiments  was  performed  using  9  mtorr  of  TFE  diluted  in 
various  pressures  of  argon  ranging  from  81  to  358  mtorr,  and  laser 
fluences  up  to  about  2  J  cm'2.  These  data  are  presented  in  Table 
I  and  Figure  4  in  the  form  of  least-squares-determined  decay  rate 
constants  and  their  associated  uncertainties.  Rate  constants 
greater  than  about  2x  103  s"1  are  unreliable,  because  of  electronic 
response  time  limitations.  At  low  fluences,  the  IRF  intensities 
are  very  low,  resulting  in  increased  uncertainty. 

For  the  higher  pressures  of  argon  and  lower  fluences,  the  gas 
temperature  after  equilibration  is  relatively  low,  but  for  high 
fluences  and  low  pressures  of  argon,  the  temperature  excursions 
are  quite  large,  as  depicted  in  Figure  5.  Thus,  the  decay  rate 
constants  may  reflect  both  dependence  on  internal  vibrational 
energy  and  dependence  on  collider  gas  temperature.  At  the  lower 
argon  partial  pressures,  collisions  between  two  excited  TFE 


(26)  Bcvington.  P  R.  ’Data  Reduction  and  Error  Analysis  for  the  Physical 
Sciences'.  McGraw-Hill.  New  York.  1969:  p  237. 


TABLE  I: 
Fla 


Infrared  FWereecence  Decay  Rates  as  a  Fa 


fluence.  J  cm' 


decay  rate,  10*  s'1 


0.21 

8 1  mtorr  of  Argon 

5.10  ±  0.65 

0.31 

5.74  ±  0.41 

0.48 

4.92  ±  0.35 

0.68 

4.68  ±  0.28 

0  86 

4.71  ±  0.22 

1.07 

4.09  ±  0.22 

1.11 

3.62  ±  0.21 

1.46 

3.67  ±0.17 

1.49 

3.25  ±  0  25 

2.00 

2.95  ±  0.18 

0.33 

1 20  mtorr  of  Argon 

8.71  ±  0.61 

0.49 

6.57±  0.49 

0.73 

6.26  ±  0  37 

0.97 

5.57  ±  0.35 

1.43 

4  65  ±  0.29 

2.25 

3.36  ±  0.2C 

030 

225  mtorr  of  Argon 

10.75  ±  1.54 

066 

10.00  ±  0.54 

0.85 

9.35  ±  0.41 

1.06 

7.69  ±  0.49 

1.33 

6  52  ±  0.31 

1  62 

6.52  ±  0.36 

1.98 

4.91  ±  0.20 

066 

358  mtorr  of  Argon 

12.60  ±  1.75 

0.92 

14.95  ±  0.98 

1.19 

11.40  ±0.81 

molecules  become  more  probable,  complicating  the  interpretation 
still  further. 

For  all  of  the  simulations,  the  functional  form  of  <A£)d  was 
assumed  to  be  the  same  as  that  used  in  previous  work: 

<A£>„  -  01 77 300|*  +  y£  (1) 

where  f),  y.  and  n  are  parameters,  T  is  the  temperature,  and  E 
is  the  vibrational  energy. 

For  computer  simulation  and  interpretation  of  the  results,  we 
arbitrarily  assumed  that  the  ( A£)d  parameters  for  collisions  with 
TFE  can  be  approximated  by  those  for  azulene.24  This  assumption 
can  be  tested  against  experiment  (see  below)  and  it  is  not  serious 
as  long  as  deactivation  by  argon  collisions  is  dominant.  The 
assumed  Lennard-Jones  collision  parameters  are  a  *  5.28  A  and 
i/K  a  323  for  TFE.  and  a  =  3.54  A  and  t/K  -  93.3  for  argon. 

For  simulation  of  the  results,  the  following  iterative  procedure 
was  used.  First,  the  data  set  for  225  mtorr  added  argon  was 
simulated,  assuming  no  dependence  of  (A £>d  on  temperature. 
Here,  the  data  set  covers  a  wide  range  of  fluences  and  the  tem¬ 
perature  excursions  are  relatively  small  due  to  the  large  amount 
of  argon  diluent.  This  simulation  gives  information  mostly  about 
the  energy  dependence  of  (A £)d,  as  distinguished  from  the  tem¬ 
perature  dependence,  because  the  initial  population  distribution 
of  vibrationally  excited  TFE  is  controlled  by  varying  the  laser 
fluence  and  the  temperature  excursions  are  not  too  great.  Next, 
the  data  set  at  81  mtorr  of  added  argon  is  simulated,  including 
the  (A£>d  temperature  dependence,  if  necessary.  Here,  the 
temperature  excursions  are  large,  providing  a  test  for  temperature 
dependence.  Last,  we  return  to  the  first  step,  if  necessary,  and 
refine  our  fitting  of  the  225-mtorr  data  based  on  a  (A £)d,  which 
depends  on  both  energy  and  temperature.  This  iterative  refinement 
procedure  results  in  an  expression  for  (A £>d  that  adequately  fits 
both  the  225-mtorr  added  argon  data  and  that  for  81-mtorT  added 
argon.  The  data  sets  obtained  at  other  pressures  of  added  argon 
provide  additional  tests  of  the  simulation. 

Following  the  prescription  given  above  for  fitting  the  experi¬ 
mental  data,  it  was  found  that  an  energy-independent  <A£)d 
simply  does  not  fit  over  the  entire  fluence  range  investigated. 
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However,  the  parameters  d  *  200  cm'1  and  energy  dependence 
y  s  0  005  gave  a  good  fit  over  the  entire  range.  In  applying  these 
same  parameters  to  the  81-mtorr  argon  data  set,  it  was  found  that 
the  resulting  agreement  was  very  close  to  the  experimental  data, 
indicating  that  <A£>( (  does  not  depend  on  temperature  in  this 
energy-transfer  system,  at  least  for  temperatures  between  400  and 
1000  K.  The  calculated  results  are  presented  in  Figure  5  as  solid 
lines  joining  the  calculated  points.  Inspection  of  the  figure  shows 
that  all  of  the  data  sets  are  adequately  represented  by  the  single, 
temperature-independent  expression: 

(l£)i  *  (200  ±  20)  +  (0.005  ±  0.002)£  (2) 

where  the  energies  are  expressed  in  cm'1  and  the  uncertainties 
are  estimated  as  the  limits  of  acceptible  simulations  (future  im¬ 
provements  in  data  treatment  may  include  better  means  for  es¬ 
timating  uncertainties).  The  simulations  seem  to  underestimate 
the  decay  rate  constants  slightly,  but  the  overall  agreement  be¬ 
tween  experiments  and  simulations  is  very  good. 

The  very  good  agreement  between  simulations  and  experiments 
indicates  that  our  starting  assumption  about  the  parameters  for 
TFE  colliding  with  TFE  is  not  seriously  in  error.  The  relative 
importance  of  such  collisions  increases  more  than  fourfold  in  going 
from  the  358-mtorr  argon  data  set  to  that  for  81-mtorr  added 
argon,  but  it  is  still  small.  The  slight  underestimate  of  the  rate 
constants  may,  however,  be  due  to  a  small  underestimation  of  the 
parameters  for  TFE-TFE  collisions. 

Equation  2  indicates  that  the  magnitude  of  <A£)d  is  signifi¬ 
cantly  smaller  than  expected  from  many  unimolecular  rate  studies 
using  other  activated  molecules.1  However,  investigation  of 
1.1,1 -TFE  by  chemical  activation27  and  by  infrared  muitiphoton 
decomposition21  studies  gave  <  A£)d  ■  350  cm'1  near  the  reaction 
threshold  (=24000  cm'1),  in  excellent  agreement  with  the  present 
results  for  1,1.2-TFE.  which  give  (A£)d  “  320  cm'1  at  the  same 
vibrational  energy.  In  the  present  study,  a  wider  variation  in 
excitation  energy  was  achieved,  permitting  a  more  direct  test  for 
the  energy  dependence  of  (A£)d. 

Comparisons  with  other  nonreactive  studies  are  not  clear-cut. 
The  results  on  deactivation  of  azulene  by  argon211  predict  <  A£}d 
=  367  cm'1  at  £  =»  24000  cm'1,  in  fair  agreement  with  the  present 
study,  but  the  energy  dependence  found  in  the  azulene  system  (7 
=  0  013)  is  much  greater  than  that  found  here.  On  the  other  hand, 
Troe  and  co-workers  have  found  no  energy  dependence  in  any 
polyatomic  energy-transfer  system  they  have  investigated,4-510  but 
the  <A£)d  values  they  found  for  argon  collider  gas  range  from 
about  230  to  420  cm"1,  in  general  agreement  with  the  present 
results. 

The  absence  of  a  temperature  dependence  is  in  good  agreement 
with  other  nonreactive  system  studies  of  energy  transfer27-10  and 
in  clear  disagreement  with  the  investigations  of  Rabinovitch  and 

I  27)  (a)  Chang,  H.  W  .  Craig,  N.  L.;  Setter.  D  W  J  Phys  Chem  1972, 
'6.  954  ib)  Marcoux.  P  J .  Siefert,  E.  E.;  Setter.  D  W.  Ini.  J  Chem  Kinet 
197*.  7.  473.  (c)  Marcoux.  P.  J.,  Setter,  D.  W  J  Phys.  Chem.  I97g,  92.  97 
idi  For  deactivation  of  other  fluoroethanea  by  gases  other  than  argon,  see 
Richmond.  G..  Setter.  D  W  J  Phys  Chem.  1990.  94.  3699 

1  28)  Jang.  J  C  .  Setter,  D  W  J  Phys  Chem  1979.  99.  3809 
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co-workers’  on  the  reactions  of  isotopically  substituted  cyclo¬ 
propanes.  The  reason  for  the  disagreement  has  not  been  identified, 
but  as  new  techniques  are  applied  to  energy-transfer  problems, 
a  better  understanding  of  the  process  will  emerge. 

It  should  be  pointed  out  that  the  heat  capacities  of  TFE  and 
argon  are  approximated  in  the  master  equation  simulations  by 
a  linear  dependence  on  temperature,  whereas  the  actual  tem¬ 
perature  dependence  is  more  complex.  Thus,  for  very  large 
temperature  excursions,  the  temperature  dependence  (or  lack  of 
it)  deduced  from  the  simulations  for  <A£)d  may  be  slightly  af¬ 
fected.  The  magnitude  of  the  error  introduced  is  related  to  the 
difference  in  the  actual  equilibrated  final  temperature  and  that 
calculated  according  to  the  linear  heat  capacity  relationship  and 
the  amount  of  heat  deposited  by  the  absorbed  laser  light.  F  - 
a  predicted  final  temperature  of  1000  K  in  a  mixture  of  9  mt, 
of  TFE  and  81  mtorr  of  A r,  the  actual  temperature  is  about  1  It*. 
K.  The  error  introduced  is  acceptably  small  at  high  temperatures 
and  can  be  neglected;  at  lower  temperatures,  the  errors  are  even 
smaller. 

In  conclusion,  this  demonstration  shows  that  the  technique  of 
infrared  multiphoton  absorption  excitation  gives  results  in  good 
agreement  with  other  methods  of  excitation,  when  the  population 
distnbutions  have  been  thoroughly  characterized.  1R-MPA  has 
the  added  advantages  over  other  excitation  methods,  however,  of 
easily  controlled  excitation  energy  (simply  by  changing  laser 
tluence)  and  of  applicability  to  many  different  classes  of  molecules. 
In  principle,  ar.y  technique  for  detection  of  excited  molecules  can 
be  combined  with  IR-MPA  for  energy-transfer  studies.  For 
example,  time-resolved  ultraviolet  absorbance  has  been  used  in 
the  past  to  monitor  the  vibrational  energy  produced  in  infrared 
multiphoton  absorption  experiments,2’  and  very  recently  it  has 
been  applied  to  collisional  energy  transfer  in  such  systems.50  In 
the  present  work,  infrared  fluorescence  has  been  shown  to  be  an 
effective  method  for  monitoring  energy  transfer. 

The  results  obtained  for  argon  collider  gas  and  excited  1,12- 
trifluoroethane  are  in  reasonable  accord  with  results  obtained  in 
other  studies  of  other  excited  molecules.  The  results  indicate  a 
small  ( A£)d  energy  dependence  and  no  detectable  dependence 
on  collider  gas  temperature. 
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A  Monte  Carlo  method  is  described  for  highly  efficient  multidimensional  integration.  The  sampling  technique  is  not  restricted 
to  hyperdimensionai  rectangles,  but  it  can  be  applied  to  more  complicated  domains.  If  known,  the  boundaries  of  an  arbitrary 
integration  region  can  be  used  to  define  the  sampling  domain,  resulting  in  Monte  Carlo  sampling  with  unit  efficiency  9s  hen 
applied  to  calculating  sums  of  states,  the  method  in  principle  can  include  all  molecular  degrees  of  freedom,  coupled  m  any 
fashion.  The  Monte  Carlo  integration  method  is  demonstrated  by  calculating  the  vibrational  sums  of  bound  states  for  HO 
and  CH.O.  based  on  their  spectroscopic  constants,  including  all  diagonal  and  off-diagonal  anharmomcmes.  The  method 
was  ioun d  to  be  highly  efficient  and  practical  for  calculations  involving  moderate-sized  molecules,  even  for  vibrational  energies 
substantially  above  the  dissociation  energy  The  resuits  show  that  inclusion  of  the  off-diagonal  anharmomcmes  produces 
a  significant  increase  in  the  calculated  sum  of  vibrational  states. 


Introduction 

Sums  and  densities  of  molecular  states  are  needed  in  calculations 
of  RRKM  ummoiecular  rate  constants,  energy-transfer  systems, 
infrared  multiphoton  absorption,  radiationless  transitions,  and  any 
system  in  which  detailed  balance  is  invoked.  Thus,  methods  are 
needed  for  calculating  the  sums  and  densities  of  molecular  states, 
regardless  of  the  nature  of  the  degrees  of  freedom  and  whether 
they  are  coupled 

In  vibrattonally  and  rotauonally  excited  molecules,  the  degrees 
of  freedom  are  strongly  coupled,  affecting  the  dynamics  of  in¬ 
tramolecular  vibrational  relaxation.  At  high  energies,  near  and 
above  reaction  thresholds,  the  coupling  can  be  very  large,  seriously 
affecting  the  sums  and  densities  of  states.  In  fact,  the  “transition 
state"— with  frequencies  different  from  the  ground-state 
molecule — is  a  device  that  is  used  to  account  for  the  effects  of 
coupling 

The  Bever-Swinehart  algorithm1  land  its  extensions  by  Stein 
and  Rabinovitch;)  revolutionized  calculations  of  molecular-state 
-.urns  iG(£))  and  densities  (p(£)  <*  dG(£)/d£).  since  its  high 
efficiency  made  exact  counts  of  states  practical  for  most  molecules. 
Despite  its  great  usefulness,  an  unfortunate  limitation  of  this 
algorithm  is  that  it  can  be  applied  only  to  separable  degrees  of 
freedom  In  principle,  the  method  of  steepest  descents  and  other 
methods  based  on  approximate  inversion  of  the  partition  function1 
can  be  used  for  nonseparable  degrees  of  freedom,  but  the  partition 
functions  are  not  known,  making  this  approach  impractical. 

Several  groups  have  applied  Monte  Carlo  techniques  to  the 
problem  of  calculating  G(E)  for  nonseparable  degrees  of  freedom. 
Monte  Carlo  techniques  are  attractive,  because  they  do  not  depend 
on  the  particular  degrees  of  freedom  or  its  coupling.*  In  a  recent 
publication.  Doll  showed  how  the  semiclassical  G(  £1  for  a  classical 
Morse  oscillator  can  be  obtained  by  application  of  “crude"  Monte 
Carlo  techniques;5  nonseparable  potentials  can  be  treated  by  the 
same  method.  Farantos.  Murrell,  and  Hajduk  applied  stratified 
Monte  Carlo  sampling  to  calculate  semiclassical  G(£Vs  and  p(£Vs 
for  coupled  vibrational  and  rotational  degrees  of  freedom.6  In 
developing  a  Monte  Carlo  version  of  transition-state  theory. 
Viswanathan.  Raff,  and  Thompson’  developing  an  interesting 
Metropolis-Markov  chain  technique  for  calculating  densities  of 
states  for  coupled,  anharmonic  potentials. 

1 1 )  Beyer.  T  .  Swmehart.  D  F  Commun  4 f  W  197J.  16.  )VS 

t3l  Stein.  S  E  .  Rabmovuch,  8  S  /  Chtm  Phis  I97J.  58.  74 J8 

i  J)  Font.  W  Theory  of  Lhtmolecular  Reactions.  Academic  New  York. 
I9'3.  p  104 

1 4)  Noid.  D  W  .  Koszykowski.  M  L  .  Tabor.  M  .  Marcus.  R  A  J  Chem 
Phil  1988,  V.  6169 

1 5)  Doll.  J  0  Chem  Phvi  Lett  t9S0.  139 

1 6 1  Faramos.  S  C  .  Murrell.  J  N  .  Hajduk.  1  C  Chem  Phis  1982.  68. 
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As  pointed  out  by  Farantos  et  al..6  much  of  the  research  in¬ 
volving  Monte  Carlo  methods  is  motivated  by  the  need  to  increase 
efficiency  The  stratified  sampling  technique  they  used  (developed 
by  Friedman  and  Wright1)  is  relatively  efficient  and  is  available 
as  a  “canned"  computer  code,  but  it  appears  to  be  cumbersome. 
Moreover,  the  “strata"  are  constrained  to  be  hyperrectangles,  even 
when  other  geometries  would  give  greater  efficiency. 

In  the  present  paper,  we  describe  a  relatively  simple  Monte 
Carlo  sampling  technique  that  gives  unbiased  estimates  and  which 
has  a  high  efficiency.  Moreover,  the  sampling  technique  can  be 
easily  applied  to  any  arbitrary  sampling  geometry.  This  feature 
may  enable  development  of  stratified  sampling  schemes  and  im¬ 
portance  sampling  schemes* 10  that  capitalize  on  the  optimum 
geometry  of  particular  problems.  Even  without  such  refinements, 
however,  the  method  is  practical  for  sums  of  states  calculations 
on  moderate-sized  molecules. 

In  testing  the  sampling  technique,  quantum  G(£)  was  calculated 
from  recent  spectroscopic  data  for  the  (nonseparable)  vibrational 
modes  of  water11  and  of  formaldehyde.15  The  results  show  that 
the  coupling  among  the  vibrations  can  significantly  increase  C(£) 
at  moderate  energies,  compared  to  the  G(£)  calculated  for  se¬ 
parable  anharmonic  vibrations.  This  result  is  particularly  im¬ 
portant  in  comparisons  of  transition-state  theory  rate  constants 
calculated  from  potential  energy  surfaces,  when  tunneling  (or  some 
other  effect)  is  invoked  to  explain  differences  from  experiments 
or  from  the  results  of  other  calculations. 

Theory 

The  sampling  algorithm  described  below  can  be  applied  to 
continuous  variables,  discrete  variables,  or  a  combination  of  the 
two.  Since  the  motivation  of  the  present  work  is  the  calculation 
of  quantum  sums  of  states,  the  descriptions  will  be  based  on 
discrete  random  variables  (quantum  numbers),  but  the  approach 
is  easily  generalized  to  any  multiple-dimensional  integral  or  sum. 

In  “crude"  Monte  Carlo*  sampling,  each  random  variable  is 
sampled  independently  and  uniformly,  resulting  in  a  uniform 
distribution  of  sample  points  throughout  the  hyperrectangle  defined 
by  the  resoective  ranges  of  the  independent  variables.  The  variable 
ranges  are  chosen  so  that  the  integration  region  lies  entirely  within 
the  sampling  hyperrectangle.  If  the  sample  points  fall  within  the 
integration  region,  they  are  counted  successes  and  the  integral 

iS)  Friedman.  J  H  .  Wright.  M  H  4CW  Trans  Sloth  Sofmare  1981. 
T.  “6 

1 9)  Hammersley.  J  M  .  Handscomb.  D  C  Monte  Carlo  Methods 
Chapman  and  Hall:  London.  1 964 

1 10)  Rubinstein,  R.  Y  Simulation  and  the  Monte  Carlo  Method,  \V  ties 
New  York.  1981 

(II)  Benedict.  W  S  .  Gailar.  N  .  Plvler.  EKJ  Chem  Phi  I  |9«6  A 

1 1 39 

( 1 7>  Reisner.  D  E  .  Field.  R  W  .  Kmsev.  J  L  Dm.  H  L  J  Chem  Ph-.  i 
1984.  SO.  5968 
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S  is  given  by  the  sampling  volume  (F,)  multiplied  by  the  ratio 
of  successes  (/i)  to  trials  (A): 

5  -  v,n/y  (i) 

If  the  integration  region  is  not  a  hyperrectangle  and  of  about  the 
same  size  as  V,.  the  efficiency  (here  defined  as  <  =  n>\)  can  be 
very  low.  resulting  in  poor  precision  of  the  Monte  Carlo  method. 

For  multidimensional  integrals,  this  poor  efficiency  can  cripple 
the  approach,  since  astronomical  numbers  of  trials  must  be  cal¬ 
culated  for  high-precision  results. 

To  improve  the  efficiency  of  the  method,  numerous  ‘variance 
reduction"  methods  have  been  developed.8"  0  The  approach 
adopted  here  is  to  choose  a  sampling  domain  that  resembles  the 
integration  region  as  closely  as  possible:  in  the  limit,  the  sampling 
domain  may  be  identical  with  the  integration  region  In  general, 
the  ideal  sampling  domain  is  not  a  hyperrectangle  of  independent 
variables  but  a  more  complicated  shape:  the  sampling  range  of 
a  given  variable  will  depend  on  values  already  selected  for  the  other 
variables.  Moreover,  the  distribution  of  sample  points  is  not 
uniform  throughout  the  sampling  domain,  necessitating  the 
calculation  of  appropriate  weighting  factors. 

For  Monte  Carlo  integration,  consider  a  multidimensional 
sampling  volume  K,  defined  in  a  space  of  discrete  variables  <e  g  . 
quantum  states).  The  gth  sample  point,  or  cell,  is  designated  by 
a  complete  set  of  variable  values  (quantum  numbers)  Equation 
1 .  for  any  type  of  sampling,  can  be  written  in  terms  of  the  total 
O  cells  and  weighted  sums  over  f[g)y,  and  St  where /(g)  and 
.V,  are  the  value  of  the  integrand  and  the  numbers  of  trials 
performed  at  the  gth  cell 

5  =  ^X>^g).vf/I>rvf  (2) 

where  the  summations  extend  over  the  C  cells.  Equation  2  can 
be  used  with  any  type  of  sampling  to  obtain  estimates  of  S 
However,  the  sampling  volume  V%  must  be  known,  as  well  as  the 
weight  function  For  complicated  sampling  domain  boundaries 
(as  in  sums  of  states  calculations).  V,  may  not  be  known. 

For  the  special  case  of  uniform  sampling,  the  sampling  domain 
is  a  hyperrectangle  of  volume  V0  containing  G  (=  oF'0)  cells  and 
the  weight  function  is  a  constant 

w0  =  l/C  =  lim  y,/y  (3) 

where  ,V  =  £.V(  is  the  total  number  of  trials  and  w0  is  the 
weighting  function 

For  uniform  sampling,  eq  2  is  an  unbiased  estimator  for  5  * 

Now  suppose  that  we  employ  weighted  sampling  with  normalized 
weight  function  ws  in  a  new  sampling  volume  V,  that  is  completely 
enclosed  within  Vy,  the  new  sampling  domain  contains  G’  cells 
(of  the  same  size  as  in  I/,),  but  we  do  not.  in  general,  know  the 
value  of  either  Vt  or  G'  The  weight  function  car.  be  expressed 
as  the  ratio  of  the  number  of  trials  ,Vf'  at  the  gth  cell  divided  by 
the  total  of  trials  <V'made  using  the  weighted  sampling  technique: 

iv.  *  lim  y.’/y  (■*) 

.v— •  * 

Using  this  identity,  we  can  connect  the  number  of  trials  at  the 
gth  cell  that  are  made  using  uniform  sampling  to  the  number  that 
are  made  using  the  weighted  sampling  method: 

lim  /lg).y,/y  =  (H/q/w ,)(/Tg).V  7 .V)  I'a) 

N  V  -* 

lim  y./y  =*  (w0/w x.v  y.v)  obt 

V \  -•  *  *  * 

For  complicated  integrands,  eq  5a  must  be  used,  but  for  simple 
determinations  of  volume,  the  function  Jig)  equals  unity  inside 
the  volume  and  zero  outside  the  volume  Thus,  m  eq  5b  we  can 
write  nt  in  place  of /fg)/Vr  where  nt  is  the  number  of  successes 
obtained  with  weighted  sampling  at  the  gth  cell  These  equations 
arc  used  to  express  the  numerator  of  eq  2  in  terms  of  nf  and  the 
denominator  m  terms  of  ,V'  =*  Kn£(,Vf7wf>.  from  eq  5 

S  «  w,nl'Z.V'~'5I/lg).Vf'/ivt  (ha) 

*  Wo F„.V  ' thbl 


where  the  summations  extend  over  the  G’  cells.  Note  that  the 

number  of  cells  in  the  hyperrectangle  is  G  *  aF0  and  Wo  »  I  G 

Thus,  w:,(/5  =  |  ’a.  and  we  obtain  an  expression  for  the  integral 
that  does  not  explicitly  depend  on  knowledge  of  U 

.9  =  (o.VVlI^g).V|ywi  (’a) 

5  =  (u.vy'EV/wr  i'b) 

Note  that  the  summations  in  eq  6  and  2  extend  only  over  the  G- 
cells  contained  in  the  new  sampling  domain,  because  there  are 
no  trials  outside  that  region.  Of  course,  not  all  cells  will  be 
sampled,  and  so  the  sum  extends  only  over  the  sampled  cells  Thus, 
we  do  not  need  to  know  the  value  of  G' 

The  precision  of  the  Monte  Carlo  estimator  is  monitored  by 
the  usual  expression*  for  the  estimated  error  of  the  mean 

</  =  i.v'uv'-  nr'Dwyrgi.vy  -sf-  <s> 

where  a  is  the  standard  deviation  of  5 

Equation  7  is  completely  general,  and  it  also  applies  to  sampling 
in  continuous  variables.  For  discrete  variables,  the  proportionality 
constant  a  is  conveniently  set  to  unity,  so  that  the  number  of  cells 
m  the  hyperrectangle  is  numerically  equal  to  I/,.  Equation  7a 
is  the  general  result,  while  eq  7b  is  specifically  appropriate  for 
calculating  volumes  in  multidimensional  space. 

Note  that  eq  7b  can  be  used  to  find  Vt,  the  volume  of  the 
weighted  sampling  domain,  if  the  boundaries  of  the  volume  are 
known  and  if  is  known.  This  fact  can  be  particularly  useful 
in  calculating  unknown  volumes  in  multidimensional  space  when 
only  the  volume  boundaries  are  known,  since  the  sampling  domain 
boundaries  can  be  set  equal  to  the  boundaries  of  the  unknown 
volume;  every  trial  will  be  a  success,  leading  to  unit  efficiency 
Nonetheless,  many  trials  must  still  be  made,  because  the  limits 
expressed  in  eq  4  and  5  are  required. 

To  utilize  eq  7.  we  need  to  define  a  normalized  distribution  of 
sample  points,  from  which  wf  can  be  calculated.  A  distribution 
of  sample  points  that  can  be  tailored  to  arbitrary  sampling  ge¬ 
ometries  is  presented  in  the  next  section 

Sample  Point  Distribution  and  Sampling  Algorithm 

The  following  algorithm  defines  a  distribution  of  sample  points 
constrained  to  an  arbitrary  sampling  domain. 

1  Choose  an  appropriate  sampling  domain  in  s-dimensional 
space  that  completely  encloses  the  integration  region,  express  it 
as 

B  =  FU-.x, . f,)  (9) 

where  eq  9  defines  the  is  -  I  )-dimensional  boundary  surface  in 
terms  of  the  variables  c.,  .  it. 

2  Starting  with  all  v,  at  their  minimum  values,  calculate  the 
maximum  value  of  c.  consistent  with  eq  9  The  range  of  r,  is  then 
given  by 

/?,  =  ('■  -  i  ;  (continuous  variables) 

R ,  =  I  +  r  -  t ;  (discrete  variables) 

3  Use  a  uniformly  distributed  random  number  L  (0  £  L  < 

1 )  to  select  a  value  of  c,  from  the  range  R , 

1 1  =  L  R,  4-  t ,  ( 10) 

(note:  truncate  r,  for  integer  variables.) 

4  Repeat  steps  2  and  3  for  each  t„  but  in  eq  9,  use  the  values 
already  selected  for  the  preceding  variables.  Continue  this  process 
until  values  for  all  variables  have  been  selected,  consistent  with 
eq  9.  This  will  result  in  the  selection  of  a  sample  point  within 
the  sampling  domain  defined  by  eq  9  Note  that  the  range  of  the 
nh  variable  is  conditional  with  respect  to  all  variable  values 
previously  selected. 

5  Calculate  the  weighting  factor.  wf,  for  sample  point  g 

w,  *  nit,-1  (in 

*■1 

*here  R,  is  the  range  of  the  #th  variable. 
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Figure  I.  Two-dimensional  sampling  solume 

6  Determine  the  value  of  f\g)  (or  whether  the  trial  is  a  success 
or  a  failure)  and  add  the  results  to  the  summations  in  eq  1 

7  Test  to  determine  whether  the  maximum  number  of  trials 
or  a  criterion  for  precision  (e  g.,  relative  standard  deviation)  has 
been  reached;  if  so.  terminate  the  calculation. 

8.  Although  not  strictly  required,  the  order  of  selection  of  the 
variables  can  be  shuffled  to  eliminate  the  variable  weighting  that 
depends  on  their  order  of  selection  (see  below). 

9  Return  to  step  I  and  initiate  a  new  trial. 

The  distribution  of  sample  points  that  results  from  this  al¬ 
gorithm  is  normalized  but  not  uniform  Furthermore,  the  variables 
are  weighted  depending  on  the  order  in  which  they  are  selected. 
This  is  most  easily  demonstrated  for  discrete  variables,  e  g  .  e, 
corresponding  to  vibrational  quantum  numbers. 

For  example,  consider  a  molecule  containing  two  separable 
vibrational  modes 
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where  E  is  the  vibrational  energy  in  excess  of  the  zero-point  energy 
and  the  w,'s  are  the  harmonic  oscillator  energy  level  spacings  This 
equation  defines  the  two-dimensional  sampling  domain  shown  in 
Figure  1  The  sum  of  vibrational  states  is  the  area  enclosed  by 
the  diagonal  line  and  the  two  axes.  (For  simplicity.  Figure  I 
depicts  the  case  when  W|  =»  but  this  condition  is  not  necessary, 
in  general )  (Note  when  anharmomcities  are  included,  the  straight 
diagonal  boundary  is  replaced  by  a  curved  line,  making  this 
problem  much  more  difficult.) 

The  sampling  algorithm  and  weighting  factor  reflect  the 
probability  of  selecting  particular  pairs  of  variables.  In  Figure 
I.  energy  E  corresponds  to  up  to  two  quanta  in  either  of  the  two 
modes  or  one  quantum  in  each  mode.  Selecting  e,  first,  the  range 
R,  =  3  (i.e.,  1 1  =  0.  1.  2);  if  the  selected  value  is  c,  =  0.  then  R: 
=  3  (i.e..  =  0.  1.  2).  but  if  c,  ■  l  or  2,  then  R:  *  2  or  1. 

respectively.  Thus,  the  probability  of  choosing  a  particular  pair 
of  values  is  ( /?, /?2)_l .  as  noted  in  Figure  1 ;  this  probability  is  just 
the  weighting  factor  wr  Note  that  the  weighting  factor  is  nor¬ 
malized: 
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As  shown  in  Figure  1.  the  distribution  of  sample  points  is  not 
uniform.  The  algorithm  tends  to  weight  the  higher  energies  more 
heavily  than  the  lower  energies.  Variables  chosen  earner  in  the 
sequence  are  more  often  selected  with  higher  quantum  numbers 
than  those  selected  later.  This  causes  no  problems  for  sufficient 
numbers  of  trials,  but  convergence  appears  to  be  faster  when  the 
selection  order  is  shuffled  between  trials.  Shuffling  eliminates 
the  weighting  due  to  selection  order,  but  the  higher  energy  points 
are  still  selected  more  often  than  the  low-energy  points.  This 
constitutes  a  fortuitous  “importance  sampling"  effect  that  is  de¬ 
sirable  in  calculating  sums  of  states,  where  the  sum  varies  more 
rapidly  at  high  energies. 

The  computational  speed  of  this  Monte  Carlo  method  depends 
on  the  efficiency  «.  which  can  equal  unity,  and  on  the  computation 


of  the  boundaries,  weights,  and  integrand  This  algorithm  cal¬ 
culates  the  weighting  factors  “on  the  fly'  and  without  knowing 
improving  the  efficiency  of  the  method,  since  weights  are 
calculated  only  for  the  points  actually  sampled  The  necessity 
of  calculating  the  boundaries  of  r,  from  eq  9  costs  some  com¬ 
putational  labor,  and  thus  simple  boundaries  are  desirable.  Other 
factors  can  also  affect  computational  speed.  For  example,  ihe 
cost  of  shuffling  the  ordering  of  the  variables  must  be  weighed 
against  the  cost  of  simply  running  more  trials.  We  have  not 
explored  all  of  these  discretionary  factors  but  have  arbitrarily 
chosen  to  shuffle  the  ordering  after  every  trial. 

Sums  of  states  for  individual  symmetry  species  are  useful  for 
various  purposes.  Several  authors  have  described  methods  for 
calculating  the  sums  of  states  for  individual  symmetry  species  in 
molecules  with  separable  degrees  of  freedom .1!~' 5  For  nonse- 
parable  degrees  of  freedom,  the  same  symmetry  species  apply.  6 
but  the  sums  of  states  will  be  affected  by  the  coupling.  The  present 
approach  lends  itself  naturally  to  calculating  sums  of  states  for 
nonseparable  degrees  of  freedom,  because  complete  sets  of 
quantum  numbers  are  selected  in  each  trial  The  symmetry  species 
of  each  vibrational  state  can  be  determined  by  application  of 
Tisza's  rules17  or  the  simplification  described  by  Lederman. 
Runnels,  and  Marcus.1’  The  number  of  successes  for  each  sym¬ 
metry  species  can  be  stored  separately  to  give  the  sums  of  states 
of  each  type. 

Sampling  Domain  for  Coupled  Vibrational  States 

For  nonseparable  anharmomc  vibrations,  the  conventional 
expression  for  the  vibrational  energy  in  excess  of  the  zero-point 
energy  is 
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where  w,  are  the  frequencies  and  Xv  are  the  anharmonicities.  This 
equation  defines  the  boundary  of  the  integration  region  for  bound 
states,  and  it  can  be  used  to  define  the  boundary  of  the  sampling 
region.  A  relatively  simple  calculated  boundary  is  obtained  if  only 
the  first  summation  term  is  retained:  this  corresponds  to  separable 
vibrations.  W  hen  the  off-diagonal  anharmonicities  are  small,  the 
“separable"  boundary  closely  approximates  the  actual  boundary 
and  it  is  a  good  choice,  but  when  the  couplings  are  larger,  the  full 
equation  can  be  used  (see  below).  To  ensure  that  the  integration 
region  is  contained  entirely  within  this  sampling  domain,  the 
boundary  is  calculated  for  the  energy  £b  =  JE.  where  d  >  I  For 
anharmonic  vibrations,  the  magnitude  of  d  will  be  of  the  order 
of  unity,  if  X„  =*  X,r 

The  most  accurate  boundary  can  be  calculated  during  the 
quantum  number  selection  process  by  taking  into  account  the 
quantum  numbers  already  selected  and  including  the  off-digonal 
anharmonicities.  This  is  done  by  using  eq  14  to  calculate  the 
energy  already  assigned  to  the  first  (k  -  I )  vibrations  and  sub¬ 
tracting  this  energy  from  £b  =  3E.  The  remaining  energy  is  then 
available  for  assignment  to  the  kth  and  subsequent  modes,  as 
follows. 

If  the  unassigned  energy  £u  is  greater  than  the  dissociation 
energy  of  the  k th  mode,  the  maximum  quantum  number  (the 
dissociation  limit)  is 


('*  mu  =  •'*  o  =  Iwt  +  2.V,t ) 
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If.  on  the  other  hand,  £„  is  less  than  the  dissociation  energy.  i\  mlI 
is  given  by 
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